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Abstract 



Anisotropic states at half-filled third and higher Landau levels are investigated in the system with a finite 
electric current. We study the response of the striped Hall state and the anisotropic charge density wave 
(ACDW) state against the injected current using the effective action. Current distributions and a current 
dependence of the total energy are determined for both states. With no injected current, the energy of 
the ACDW state is lower than that of the striped Hall state. We find that the energy of the ACDW state 
increases faster than that of the striped Hall state as the injected current increases. Hence, the striped 
Hall state becomes the lower energy state when the current exceeds the critical value. The critical value 
is estimated at about 0.04 - 0.05 nA, which is much smaller than the current used in the experiments. 
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Chapter 1 

Introduction 



Since the discovery of the integer quantum Hall effect by von Klitzing et al. [l| in 1980, a two-dimensional 
(2D) electron system in a perpendicular magnetic field has attracted much attention for more than 
a quarter of a century. In 1980, von Klizing et al. have measured the Hall conductivity of the 2D 
electron system in strong perpendicular magnetic fields at low temperatures, and have found that the 
Hall conductivity exhibits plateaus just at integer multiples of e 2 /h as the magnitude of the magnetic field 
is varied, where e(> 0) is the electron charge and h is the Plank constant. This is the integer quantum 
Hall effect (IQHE), and after its discovery, a 2D electron system subjected to a perpendicular magnetic 
field is called a quantum Hall system. The integer quantization of the Hall conductivity is closely related 
to a Landau level quantization of an energy spectrum of free electrons as follows. In classical mechanics, 
a free electron performs a cyclotron motion in a uniform magnetic field B. Since this cyclotron motion 
corresponds to a harmonic oscillation with the cyclotron frequency lo c = eB/m e (m e : electron mass) 
in quantum mechanics, the energy spectrum is quantized to Ei = fku c (l + 1/2) (I = 0, 1,2, . . .). This 
quantized energy level is called Landau level (LL) named after Landau who solved this problem for the 
first time. The number of states in each LL is proportional to B, and when the magnetic field is strong 
and the lower few LLs are fully occupied by electrons, the Hall conductivity is quantized. Thus, the IQHE 
is closely related to the LL quantization. In order to explain the emergency of the plateau structure, 
a localization effect of electrons trapped by impurities should be taken into account and this provides 
another interesting topic. For more details about the IQHE, sec the reviews by Ando et al. Q and 
Stormer 

Only two years later after the discovery of the IQHE, another surprising phenomenon, a fractional 
quantum Hall effect (FQHE), was found. In 1982, working with much higher mobility sample, Tsui et 
al. [H have discovered the fractional quantization of the Hall conductivity. The fractional quantizations 
have been observed around characteristic partial fillings in the lowest LL. Around these fillings, electron- 
electron interaction plays a predominant role, which results in a unique collective phenomenon. Among 
a number of theoretical studies done to solve this problem, it was Laughlin [B| who has succeeded in 
discovering a Laughlin's variational wave function which explains the fractional quantization at the char- 
acteristic fillings. The fact that the Laughlin wave function is a good variational wave function implies 
that the ground state at these fillings is a translationally invariant liquid state. Now, various studies 
support this conclusion. After the discovery of the Laughlin's wave function, his novel work has evolved 
into a concept of a composite Fermion which is an electron bound to an integer number of magnetic flux 
quanta. Composite Fermion is based on the observation by Jain Q that the Laughlin's wave function 
could be viewed as integer quantum Hall states of composite Fermions with an even number of flux 
quanta, in a very natural way. On the basis of this prominent idea, most experimental results for the 
FQHE can be explained naturally, which supports the idea of the composite particles experimentally. 
The composite Fermion has provided exciting fields of theoretical studies since its discovery by Jain. For 
more details, see the reviews by Jain and Kamilla Q . 
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CHAPTER 1. INTRODUCTION 



There are many interesting topics other than the IQHE and the FQHE, especially in higher LLs 
where various electron solid phases exist. In 1996, Koulakov et al. [|[ have studied a possibility of charge 
density wave (CDW) states in higher LLs, which include bubble states (CDW states with two or more 
electrons per unit cell) and striped Hall states (unidirectional charge density wave states), within the 
Hartree-Fock (HF) theory. They have found that CDW and bubble states are energetically favorable in 
higher LLs, and in particular around half fillings, striped Hall states are energetically favorable. Indeed, 
in 1999, working with ultra-high mobility samples at low temperatures, Lilly et al. Q and Du et al. 
[To| | navc reported the highly anisotropy in the longitudinal resistivity, which is interpreted as due to a 
formation of stripe states or anisotropic charge density wave (ACDW) states. Away from half-fillings, 
more precisely around v = 4 + 1/4, 4 + 3/4, reentrant integral quantizations of the Hall conductivity 
have been observed 13, U|- This reentrant integer quantum Hall effect (RIQHE) is interpreted as due 



to a collective insulating property of a bubble state pinned by impurities. Although these intriguing 
phenomena in higher LLs have been studied intensively over ten years, many issues still remain to be 
solved. Among them, we focus on the highly anisotropic states. While experimental features of the 
anisotropic states suggest that the anisotropic states are the striped Hall states, the ACDW state has 
a slightly lower energy than the striped Hall state within the HF theory with no electric current. This 
has been an enigma as to the anisotropic states. In experiments of the anisotropic states, current is 
injected. This effect has not been taken into account in the previous calculations. We focus on this fact 
and investigate the effect of the injected current on the striped Hall states and the ACDW states by 
calculating response functions against the injected current [12j. This is one of main topics of this thesis. 
The details are given in Chapter 0] 



Another main topic of this thesis is a calculational method for periodic states in the quantum Hall 
system by means of a von Neumann lattice (vNL) basis 13|, 14 1. Most theoretical studies of the quantum 



Hall system start by expansion of the Hamiltonian using a complete set of eigenstates of free electrons in 
a magnetic field. There are mainly three different complete sets to be used according to the purpose: 

1. Eigenstates of a magnetic angular momentum operator 

2. Eigenstates of a magnetic translation operator in one-direction 

3. A complete set of a vNL basis which we use aggressively in this thesis 

All of them are closely related to a magnetic flux 0o = h/e, and </>o is closely related to a magnetic 
translation group as follows. As has been mentioned, eigenenergies of electrons in a uniform magnetic 
field B are quantized to LLs and each LL has a large number of degeneracy given by A^ = SB/4>q, 
where S is a total area of the 2D system. This number of degeneracy is equal to the number of magnetic 
flux quanta penetrating the system, which means that one state in each LL occupies an area penetrated 
by one magnetic flux. This fact can be viewed from a point of view of a magnetic translation. In the 
case of no magnetic field, translation operators commute with each other and they commute with a free 
Hamiltonian, too. Therefore, a momentum is well-defined as a good quantum number. On the other 
hand, in the presence of a uniform magnetic field, previous translation operators are not commutative 
with the free Hamiltonian anymore. This is because the vector potential generating the magnetic field 
is not uniform in coordinate space even if a magnetic field is uniform. To obtain the proper translation 
operator, the previous translation operator should be modified to a magnetic translation operator which 
is commutative with the free Hamiltonian in the uniform magnetic field. After we obtain the proper 
translation operators, two linearly independent magnetic translation operators are still non-commutative 
with each other for an arbitrary translation. Only discrete translation operators are commutative with 
each other, which span a magnetic translation group [la |. For these discrete operators, a momentum is 
well-defined as a good quantum number, and the momentum is reduced to a magnetic Brillouin zone 
owing to the discrete translation symmetry. Actually, the vNL basis is an eigenstate of these discrete 
magnetic translation operators so that it has a momentum as a good quantum number and periodicities 
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in two linearly independent directions. The area spanned by the two different unit magnetic translation 
operators is given by S a = 4>o/B = h/eB which corresponds to an area penetrated by <po. Hence, the 
area occupied by one electron and a magnetic translation group are closely related to each other. 
Three complete sets mentioned above are discretized as one state occupies an area S a as follows: 

1. For the complete set of cigenstatcs of a magnetic angular momentum operator, the eigenstate is 
characterized by the eigenvalue of a magnetic angular momentum operator h(l — to), where I is 
the LL index and m is a nonnegative integer. The probability density of the eigenstate with to is 
localized on the circumference with a radius r m = y^m/ira with a = ^Jh/eB, which gives S a per 
eigenstate. 

2. For the complete set of eigenstates of a translation operator in one direction (which is set to y), while 
the probability density in the y-dircction is uniform, the density in the x-direction is localized at 
discrete coordinates with the interval Ax = <Pq/BL v , where L y is the system size in the y-direction, 
which gives S a per eigenstate, too. 

3. For the complete set of the vNL basis, the probability density is localized at lattice sites with an 
area of the unit cell given by S a , which naturally gives S a per eigenstate. 

Thus, all of three basis sets are closely related to <fio which is closely related to the magnetic translation 
group. As for the vNL basis, the periodicity of the vNL can be deformed as long as the area of the unit 
cell is kept to S a . From this point of view, we develop the vNL formalism to analyze periodic states, 
which include striped Hall states, ACDW states, and Bubble states. The vNL basis is the most suitable 
basis among three sets of basis mentioned above since it can be deformed as the periodicity of the vNL 
coincides with the periodicity of the periodic state. The details are given in Chapter [3] The study of the 
highly anisotropic states in this thesis is one of applications of this formalism. There will be some other 
applications of this formalism to study periodic states in the quantum Hall system, which is one of our 
future plans. 

This thesis is organized as follows. In Chapter [21 a von Neumann lattice formalism is reviewed. Some 
useful relations and features peculiar to the von Neumann lattice formalism are also reviewed. In Chapter 
[31 a calculational method by means of the von Neumann lattice formalism is developed to study periodic 
states in the quantum Hall system. In Chapter [H effects of an injected current on highly anisotropic 
states are studied using the von Neumann lattice formalism developed in Chapter [31 Some concrete 
calculations are presented in Appendixes. 



CHAPTER 1. INTRODUCTION 



Chapter 2 



Review of the von Neumann lattice 
formalism 

In this chapter, we review the von Neumann lattice formalism which is a powerful tool to study periodic 
states in the quantum Hall system as has been mentioned in the previous chapter. The von Neumann 
lattice basis is an eigenstate of commutative magnetic translation operators, which spans a Hilbert space 
in the quantum Hall system in conjunction with a Landau level space. The von Neumann lattice basis is 
first introduced as a complete set of coherent states. Since they are coherent states, two different states 
are not orthogonalizcd. However, by Fourier transforming and normalizing them, we obtain a discrete set 
of orthonormal von Neumann lattice basis as a complete set of one-particle eigenstates. We review this 
procedure and show that the obtained basis is indeed the eigenstate of commutative magnetic translation 
operators. Some useful relations and features peculiar to the von Neumann lattice basis are also reviewed. 

Note that in this thesis, we only consider the rectangular von Neumann lattice basis for simplicity. 
However, the oblique von Neumann lattice can be obtained by a similar procedure. For details, see 
Ref. 0. 

2.1 Orthonormal von Neumann lattice basis 
2.1.1 Landau level quantization 

Let us consider a two-dimensional (2D) electron system in the x-y plane subjected to a perpendicular 
magnetic field B = (0,0, B). In the absence of electron-electron interactions, electrons just perform a 
cyclotron motion with a cyclotron frequency ui c = eB/m e . Since this cyclotron motion corresponds to 
a harmonic oscillator in quantum mechanics, one expects that the energy level is quantized to Ei = 
huj c (l + 1/2) with I = 0, 1,2, . . . . We confirm this quantization by a microscopic analysis below. 

The cyclotron motion is represented by a guiding center coordinates X = (X, Y) and the relative 
coordinates £ = (£, 77) as 

x = X + i, y = Y + r,, (2.1) 

^ = ^-(-ihdy + eAy), r] = -^-(-ihd x + eA x ). (2.2) 

where B = V x A(x) and e > is the electron charge. Here, X and £ satisfy the following commutation 
relations, 

[X, Y] = -[£, r,} = i A [X, £] = [X, r,] = [Y, = [Y, V ] = 0. (2.3) 



r 
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By using these coordinates, the one-particle free Hamiltonian Hq = (—ihV + eA(x)) 2 /2m e is expressed 
by 

Ho = ^(e + V 2 )- (2.4) 
Next, we introduce new operators La, L a , Lb, L b as 



1 ,„ 1 



L B = -^-(X + iY), L B = ^-(X-iY), (2.5) 



where Ib = y h/eB is the magnetic length. Since these operators satisfy the following commutation 
relations, 

[L A , L\] = [L B ,L B ] = 1, [La, L b ] = [L\, L b ] = [L A , L g ] = [L A , L B ] = 0. (2.6) 

new operators can be regarded as independent two kinds of creation and annihilation operators. Using 
these operators, we rewrite the one-particle free Hamiltonian Hq as 

H Q = hw c (L A L A + \). (2.7) 

This is a Hamiltonian of a harmonic oscillator with a frequency lu c . The eigenstate and eigenvalue of Hq 
are given by 

|/,) = ^|/ >, H \f l )=E l \f l ), (2.8) 

where |/o) satisfies LaI/q) = 0. Hence, the naive expectation with respect to the quantization of the 
energy level is confirmed. The energy level labeled by I, Ei = hu> c {l + 1/2), is called the Ith Landau level 
(LL). 



2.1.2 Eigenstate of the magnetic angular momentum operator 

In the above discussion, the guiding center coordinates are not used to derive the LL quantization. Since 
the free Hamiltonian Hq is commutative with the guiding center coordinates, each LL has degeneracy 
coming from the eigenstate related to the guiding center coordinates. The most simple one is the eigen- 
state of the operator L b Lb- This eigenstate is given by 

(j) )M 

\J M ) = ±-?J=-\J ), L B L B \J M ) = M\J M ) (M = 0,1,2,...), (2.9) 

where Lb\Jq) = 0. The direct product \fi) <E> \ Jm) = is a simultaneous eigenstate of the free 

Hamiltonian Hq and the operator J = (e£?/2)(£ 2 + r/ 2 — X 2 — Y 2 ) = H(L a La — L b Lb), and spans 
a Hilbert space of a one-particle state. Actually, the operator J is a magnetic angular momentum 
operator since it commutes with the free Hamiltonian, and when the magnetic field is set to zero, it takes 
J(B = 0) = —ih(xdy — yd x ) which is an angular momentum operator in the absence of a magnetic field. 
Here, the fact that A(x) is proportional to B is used. The electron density of the eigenstate with M 
is localized on the circumference with a radius rM = M/ir a with a = y^h/eB, which gives the area 
S a = h/eB = 4>o/B (4>q — h/e is a magnetic flux quantum) per eigenstate as has been mentioned in the 
previous chapter (see also Appendix lA.ip . Note that each LL has the degeneracy factor N$ = BS/<pQ (S 
is a total area of the 2D system) , which means that the number of degeneracy in each LL is equal to the 
number of magnetic flux quanta penetrating the 2D system. 
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2.1.3 Coherent state 

There is another eigenstate constructed by the guiding center coordinate, that is, a coherent states of the 
guiding center coordinates, 

TL 

(X + iY)\a mn ) = z mn \a m , n ), z mn = a{r s m + i—), (2-10) 

r s 

where m and n are integers. In coordinate space, these coherent states are localized at the rectangular 
lattice point a(mr s ,n/r s ), where a = ^Jh/eB is a lattice constant, and r s is an asymmetry parameter 
of the unit cell (see Appendix lA.lj) . The magnitude of a is of the order of tens of nano meter for a few 
Tesla of magnetic field. Since the number of these coherent states in each Landau level is equal to N^, 
the direct product \ fi) £g> |a mi „) = |^,N) with N = (m,n) forms the complete set of a one-particle state. 
Note that the completeness of this coherent set is ensured mathematically 0, | 17l Il8j and this set is 
called von Neumann lattice (vNL) basis [H, H3l- 

Since [Lb,L b ] = 1, we can take Lb = d/dL B . Thus, 

LsWm.n) = T-\a m ,n) = \0L m .n), ( 2 -ll) 

oV B v2/s 

and \a m ,n) is given by 

\a m ,n) = C m , n e^^ l ^ L B |a 0)0 ), (2.12) 

where C m . n is a normalization constant. If we set {a m>n \a m>n ) = 1, then |C m .„| = e - ^' 2 " 1 "' l 2a with 
an arbitrary phase factor. We use e i7r (™+"+ m ' 1 ) as a phase factor of C rn . n throughout this thesis and set 
a = 1 unless otherwise stated. 

The eigenstate |a„ i;rl ) is not an orthogonal basis since it is a coherent state. However, the inner 
product, 



gi7r{ (m— m') + (n— n')+(m— m')(n— ( 7r /2){(^i— m') 2 r^ + (n— n' ) 2 / r 2 } ^ ^^) 



is a function of the difference, (m — m') and (n — n'), so that the Fourier representation of |a m;ra ) becomes 
an orthogonal basis. Indeed, by using the momentum representation of \a m ,n), 



y p ) =J2 emPx+mPv \ a m,n), (2.14) 
the inner product of \a p ) is calculated as (see Appendix IA.2|) 



P |<V) - /3*(p)/3(p')(27r) 2 J2$ 2 (P - P' - 27rN)e i *(P' N ) > 

irl ) , (2.15) 



N 

„2 



/3(p) = (V^r^e^^^e, Px Y sPv 

where 0i(v|r) = i En=-oo e l7r,l+l7rr ( ,l - 1 / 2 ) 2+i7rl '( 2 "^ 1 ) is a Jacobi's theta function of the first kind, cf>(p, N) = 
7r(m + ra) — p x ri, and we set a = 1. Thus, we obtain the orthonormal vNL basis by normalizing |a p ) as 

IA>) = ^)l«p>- ( 2 - 16 ) 

From the property of the theta function, /3(p) and \ f3 p ) obeys the following nontrivial boundary condition, 
/3(p + 2ttN) = e^ N >/3(p), |/3 p+27rN ) = e-^P^|/3 p ). (2.17) 
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Using this boundary condition, we obtain the momentum space reduced to the Brillouin zone (BZ), 
— 7r < p x ,Py < ir, hence, the Hilbert space of a one-particle state is spanned by the state \l, p) = |/z)<£>|/3 p ). 
Owing to the boundary condition (|2.17[) . |/3 P ) satisfies 

(/? P |/V) = (2tt) 2 Y,£ 2 (P-P'- 27rN) e ^P< N ). (2.18) 

N 

The left hand side and the right hand side of Eq. (|2.18|) obey the same boundary condition. If we restrict 
the momentum range to the first Brillouin zone, then Eq. ([2~Tgjl becomes (/3 P |/V) = (2tt) 2 ^(p-p') which 
expresses a usual orthonormal relation. 

2.2 Magnetic translation group and von Neumann Lattice basis 
2.2.1 Magnetic translation group 

The translation operator in the system with zero magnetic field is given by 

T(6) = e lS - p/h (2.19) 

where p = —ihV. Indeed, if T(8) acts on a one-particle wave function tp{x), then tp(x.) is translated 
to ^(x + 8). The translation operator T(d) commutes with the one-particle free Hamiltonian Hq(B = 
0) = p 2 /2m e , and T(6 X ) and T(8 y ) are commutative for an arbitrary S x and S y , where S x = S x e x , 
8 y = 5 y e y , and e x and e y are unit vectors in the x and y directions, respectively. Hence, we can take 
{Hq, T(S x ) 1 T(S y )} as a set which can be diagonalized simultaneously. The situation becomes different 
when a magnetic field is applied. Let us consider the 2D electron system in a uniform magnetic field 
B = (0,0, B). In this case, the one-particle Hamiltonian is given by Ho = (p + eA(x)) 2 /2m e and the 
electron performs a cyclotron motion. Clearly, T(8) does not commute with Hq since the vector potential 
A(x) depends on x even if the magnetic field is uniform. In order to obtain the proper translation 
operator, magnetic translation operator, which commutes with Hq, the gauge transformation which brings 
A(x + £) back to A(x) is required in conjunction with T(S). This magnetic translation operator is given 
by 

T(S) = e lS - K/h , K = p + eA(x) - eB x x. (2.20) 

Since K is rewritten as K = eB(Y, — X) and Hq is rewritten as Eq. (|2.4|) . it is clear that T(S) indeed 
commutes with Hq. 

The magnetic translation operator T(S) includes the spatial translation and the gauge transformation. 
To make this point clear, we divide T(8) into two parts. Using the Cambcll-Hausdorff formula, 

cxp A cxp B = cxp{A + B + ^ [A, B] + [A - B, [A, B]] + . . . ) , (2.21) 

we rewrite T(S) as 

T(S) = e ie / 2 e iS - e ( A -' Bxx )/ h e iS -p/ h j (2.22) 

where iO = [id ■ e(A(x) — B x x)/S, id ■ p/H] is constant and the vector potential A(x) is assumed to be a 
linear function of x. This expression clearly shows that T(S) consists of the spatial translation operator 
T(S) and the gauge transformation operator e lS / 2 e tS - e ( A ~' Bxx ) / h _ 

The magnetic translation operators do not commute with each other for arbitrary two translations 
owing to the non-commutativity of K x and K y , i.e., [K x , K y ] = iheB. For simplicity, we consider two 
linearly independent magnetic translation operators, T(S X ) = e lSxK ' c ^ h and T{S y ) = e lSyKy / h . Using the 
commutation relation of K x and K y , we obtain 

f{S x )f{S y ) = f{8 y )f{8 x )e- ieB5 * 5 *> h . (2.23) 



2.3. VON NEUMANN LATTICE FORMALISM 
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Clearly, T(S X ) does not commute with T(S y ) for arbitrary S x and 5 y . However, if $ = B8 x 5 y is equal 
to 2mrh/e with an integer n, T(S X ) commutes with T(d y ). Here, notice that <!> represents the magnetic 
flux penetrating the rectangular area which is spanned by 5 X and S y , and </>o = 2?:h/e = h/e represents 
a magnetic flux quantum. This implies that T(d x ) and T(d y ) are commutative only when the magnetic 
flux penetrating the area spanned by 6 X and 6 y is equal to integral multiples of a magnetic flux quantum 
4>q . This statement can be easily generalized for arbitrary two linearly-independent magnetic translation 
operators. Hence, in the system with a uniform magnetic field, we can take {Hq, T(8\), T(<52)} with <5i 
and 82 which satisfy the relation B • (8i x £2) = >fion, as the simultaneously diagonalizable set. T(8i) and 
T(<$2) form a group called magnetic translation group [lij ]. 



2.2.2 Eigenstate of the commutative magnetic translation operators 

The vNL basis is a simultaneous eigenstate of the set {Hq, T(r s a, 0), T(0, a/r s )}, so that it is an irreducible 
representation of magnetic translation group. Let us check this statement in what follows. 

First, T(r s a,0) and T(0,a/r s ) are the magnetic translation operators along the x and y direction of 
a rectangular unit cell, respectively. Since the magnetic flux penetrating the unit cell is equal to (j>o, 
T(r s a,0) and T(0,a/r s ) are commutative. Next, when T(r s a,0) and T(0,a/r s ) act on |a m! „), \a m _ n ) is 
translated to the nearest neighbor states, |a rn _i !n ) and |a mi „_i) up to a phase factor, respectively, that 
is, 

T(r s a, 0)|a m n ) = -|a m _i „), T(0, — )\a m n ) = -|a m n _i). (2.24) 

r s 

Equation (|2.24|) gives the action of T(r s a, 0) and T(0,a/r s ) on |/3 P ) by 

f(r s a,O)\0 p ) = -e^|/3 p ), T(0, — )|/3 P ) = -e^|/3 p ), (2.25) 

rs 

hence, the vNL basis |/3 P ) is a simultaneous eigenstate with eigenvalues — e lPx and — e Wy for T(r s a, 0) 
and T(0, a/r s ), respectively. Note that the minus sign in Eq. (|2.25|) is a convention and we can take a 
different phase factor instead of the minus sign in Eq. (|2.25p . in which case the vNL basis \(3 P ) has a 
different normalization factor lj|, 2(| . 



2.3 von Neumann lattice formalism 

In this section, we derive some useful relations in the quantum Hall system in the vNL formalism. We use 
the natural unit (h = c = 1), set a = 1, and ignore the spin degree of freedom, unless otherwise stated. 

2.3.1 Hamiltonian of the quantum Hall system 

In the second quantized form, the Hamiltonian of the quantum Hall system is given by 

H = JC + V, (2.26) 

K = I rfW(x) ( ~' V + eA W)2 W(x), V= - / d 2 xd 2 x' : p(x)V(x - x'Wx') :, (2.27) 
J 2m e 2 J 

where ^(x) is an electron field operator, p(x) = ^(x)'I'(x) is a density operator, colons represent a 
normal ordering with respect to creation and annihilation operators, V(x) = 47re 2 /e|x| is a Coulomb 
potential, and e is the dielectric constant. The electron field operator is expanded by the vNL basis as 

*( x ) = E / BZ (2^ Mp)(x|/ ' p> ' (2 - 28) 



1=0 
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where bi (p) obeys the boundary condition, 

Mp + 2ttN) = e^P' N )6 ; (p), (2.29) 
and satisfies the following anti-commutation relation, 

{6,(p), bl(p')} = ^(2tt) 2 5 2 ( P - p' - 27rN) e ^(P' N ). (2.30) 

N 

The Fourier transform of the density operator p(k) = J d 2 x ^T( x )^ r (x)e is written as ( Appendix IA.3[) 

<°( k )=E / 7^^(p)6«'(p-k)/{J,,(k)e-( fr -/**)^^-*»/'--), (2.31) 
i,v J BZ * 71 > 

where / ; ° ; /(k) = (fi\e (the explicit form is given in Appendix IA.4|) and k = (r s k x ,k y /r s ). Note 

that the integrand of p(k) is invariant under the transformation p — > p + 27rN. Substituting Eqs. (|2.28l) 
and (I2.3ip into Eq. (|2.26p , we obtain the kinetic term K. given by 

/C = E i? ' j ( z (0^(P^(P) J E l= u c {l + \l (2.32) 

and the Coulomb interaction term given by 

V = - J ^ : p(k)V(k)p(-k) :, ^(k) = -j^- (k^O), 9 2 = — , (2.33) 

where V(0) = owing to the charge neutrality of the system. Using the new density operator defined by 

«,r(k)=/ ^bl(p)b,(p-k)e-^/^ k ^- k ^\ (2.34) 
Jbz ( 27r ) 

we write p(k) as /?(k) = / ; ° ; , (k)^// (k) and obtain 

V = \\j^k E ^( k )/°^( k )/° 3 ,u(-k):^^(k)p i3>i 4(-k):. (2.35) 



2.3.2 Hartree-Fock Hamiltonian 

The Hamiltonian of the quantum Hall system takes a simple form in the Hartree-Fock (HF) approximation 
by virtue of the magnetic field. Both the Hartree term and the Fock term are proportional to the density 
operator pij/fe) (Appendix IA.5[) . The Coulomb interaction term V is approximated by Vhf — (Vhf)/2, 
where Vhf is given by 

V hf= E J ^ v Z 2 ,hM^)(Ph^))PlsjM (2-36) 

with the HF potential 

<L, i3 ,* 4 (k) = U(k)/°, 2 (-k)/°, 4 (k) - J ^^(kO^^J-kO^^CkOe-™^^-^^), (2.37) 

where k = (k x /r s ,r s k y ). In Eq. (|2.37jl . the first term and the second term in the right hand side represent 
the Hartree term and the Fock term, respectively. 



2.3. VON NEUMANN LATTICE FORMALISM 



13 



If the Hamiltonian is projected to the Ith LL, the kinetic term is quenched and the Hamiltonian is 
given only by the Coulomb interaction term projected to the Ith LL. In this case, the Ith LL projected 
Hamiltonian is given by 

n(l) =\l (£H k) : «( k )a(- k ) : • ( 2 - 38 ) 

with 

v l (k) = V(k)[F l (k)} 2 , F,(k) = / I » (2.39) 
pi(k)=/ -^6|(p)6 i (p-k)e^ (lrs/4 " ) ^ (2p! '^ /rs) . (2.40) 

We call pik) a projected density operator. 

The projected HF Hamiltonian is given by H^ F = Tl^p — (Ti. l ^ F )/2, where 

n£ F = J ^vrfafai-typifc), (2.41) 

with 

v? F (k) = u,(k) - / -^^(k^eW 2 *^**-^*-). (2.42) 
7 (2tt) 2 

These notations will be used in the following chapters. 



2.3.3 Density operator 

The projected density operator (|2.40[) is noncommutative. Indeed, the commutation relation of projected 
density operators is calculated using Eq. (|2.30|) as 

fa(k),Pl0O] = -2isin( (k * k ' } ' ) ptik + k'). (2.43) 

Owing to this noncommutativity, various phases are realized in the quantum Hall system. 

The operator of the total number of electrons is given by iVtotai = / d 2 xp(x) = p(k = 0). Substituting 
Eq. (|2.31|) into this expression, we obtain 

#totai = W -p^b\(p)bi,(p){f l \f l ,)=J2 I tSU(p)&j(p)- (2-44) 

j v JBZ V Z7r ) j J BZ (Z7TJ 

Hence, in the LL projected space, the density operator in the momentum space is given by bj(p)bi(p). 



2.3.4 Magnetic field in momentum space 

In the vNL formalism, nontrivial phase factors appear in expressions such as the commutation relation 
of the field operator Eq. (|2.30|) and the density operator Eq. (|2.31|l . This nontrivial phase factor can be 
interpreted as due to the magnetic field in momentum space. On the vNL, the momentum is defined in 
the Brillouin zone and has a periodicity 2n in both p x and p y directions. Hence, the momentum is defined 
on a 2D torus. Since the field operator 6j(p) obeys a boundary condition 6((p + 27rN) = e l< ^ p ' N '&;(p), for 
one period, the phases e %lx in the p x -direction and e l7r_lPx in the p H -direction are obtained. If these phases 
are interpreted as due to the Aharonov-Bohm phase caused by the magnetic field in momentum space, 
the phase e l7r can be regarded as due to two magnetic fluxes with the magnitude 27r shown in Fig. 12. H and 
the phase e~ lVl can be regarded as due to the magnetic field perpendicular to the surface of the torus 
with the magnitude —l/2ir. Indeed, the expression of the density operator Eq. (|2.3ip is rewritten as 
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^271 




Figure 2.1: Two magnetic fluxes in the 2D torus in momentum space. 

=E J BZ ■§^ b ^ h '^ %fl,l'to CX P k A (P') ' d P'| . ( 2 - 45 ) 

where the p' integral is a line integral on the straight line from p to p — k and A(p) = — (a 2 /27r)(— p y , 0, 0) 
which corresponds to a vector potential in momentum in Landau gauge. The magnetic field obtained 
from this vector potential is B = (0, 0, — 1/2-k). The phase factor in Eq. (|2.45[) is a gauge connection 
form p to p — k in order to keep the integrand of p(k) invariant against the gauge transformation in 
momentum space 

6,(p) -> e lHp) h(p), A(p)-A(p)-VA(p). (2.46) 

Hence, the nontrivial phase factor in the expression of the density operator Eq. (|2.31|) can be interpreted 
as due to the magnetic field in momentum space consistently. While the vector potential takes the form 
in Landau gauge, other gauges can be taken by the gauge transformation in momentum space. 



Chapter 3 

Periodic states in the quantum Hall 

system 



In this chapter, we develop a formalism based on the von Neumann lattice (vNL) basis to treat periodic 
states in the quantum Hall system, which include charge density wave (CDW) states or Wigner crystal 
states, bubble states, and striped Hall states. Wigner crystal states or CDW states have been observed 
at low partial fillings and e.g., the Landau level filling factor v = 1/5 in the lowest Landau level (LL) 



2l| , whereas bubble states and anisotropic charge density wave ( ACDW) or striped Hall states have been 



observed at, e.g., partial filling factor v* ~ 1/4, 3/4 in higher LLs for the bubble state [lCJ,[ll| and v* ~ 
1/2 in the higher LLs for the ACDW or striped Hall state 0, [l(J. The vNL basis is the most suitable 
basis to study these periodic states consistently We first review the background of this issue, then we 
develop the vNL formalism for periodic states in the quantum Hall system. 



3.1 Background 

In the quantum Hall system, energy levels split into LLs owing to the cyclotron motion of electrons and 
the energy difference between the nearest LLs is given by hu> c . If the magnetic field is so strong that the 
energy difference hw c is much larger than the typical order of the Coulomb interaction e 2 /Airels, the LL 
mixing effects can be neglected. In this case, it is enough to consider only the system projected to the 
uppermost partially-filled LL where the kinetic term is quenched and the Hamiltonian is given by only the 
projected Coulomb interaction term. This projected Coulomb interaction term is totally different from 
the unprojected Coulomb interaction term in that the projected density operators are noncommutativc 
with each other as seen in Eq. (|2.43p . If we neglect this noncommutativity, an expected ground state in 
this system is a Wigner crystal state. A Wigner crystal state is a state where electrons form a triangular 
lattice through the Coulomb interaction so that the Coulomb energy becomes minimum. When the partial 
filling factor is small enough, the triangular Wigner crystal is indeed the ground state. This is because 
for low partial fillings, the superposition of wave functions is negligible and the problem can be treated 
as a classical one. The situation totally changes when the superposition of wave functions becomes large 
and the problem cannot be treated as a classical one anymore. In this case, if we exclude a possibility of 
liquid states or fractional quantum Hall states, naively expected ground states are "CDW siaies". The 
CDW states are obtained as a self-consistent solution within the Hartree-Fock (HF) approximation by 



proper treatment of the noncommutativity of the density operators [24|. In the low partial filling limit, 
the self-consistcntly obtained CDW state coincides with the triangular Wigner crystal state. However, 
as the partial filling factor approaches the half-filling, bubble states which are CDW states with two 
or more electrons per unit cell, ACDW states, and striped Hall states which are unidirectional charge 
density wave states, become ground states. To treat these states, some theoretical approaches have been 
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developed in early studies. 



3.1.1 Early studies 

Among early studies of CDW states in the quantum Hall system 22, 23|, 24], 25 1, the self-consistent HF 
solution was first studied by Yoshioka and Fukuyama [3]. They focused on the triangular and square 
CDW state in the lowest LL and calculated the energy dispersion and the total energy of the ground state, 
neglecting the higher harmonics of the CDW. The calculation was done by means of the Green function. 
In their calculation, first, the Hamiltonian is expanded by the eigenstate of the magnetic translation 
operator in one direction, and then, the Green function was ingeniously introduced to solve the problem 
self-consistently. Although their formalism is a bit complicated, the self-consistent calculation including 
the noncommutativity of the projected density operator is performed properly. Later, Yoshioka and Lee 
[26| developed a simpler formalism and calculated the same problem including higher harmonics of the 
CDW, in a different way. Using the same basis as the previous one, they more directly obtained the 
solution by diagonalizing the HF Hamiltonian self-consistently. Their results are essentially the same as 
the previous one except for small corrections. 

There is a totally different approach to obtain the same results first developed by Cote and MacDonald 
[27I ] , which is an approach using an equation of motion. They calculated the total energy and the mean 
value of the density operator for the CDW ground state by solving the HF equation of motion numerically, 
and obtained the same results as the previous one. As pointed out in their paper (27j . it is not possible to 
obtain the one-particle energy dispersion with this approach in contrast to the previous studies [24l . |26| . 
However, their main topic was to investigate collective modes of the CDW state, and for this purpose, 
it was enough to know the mean value of the density operator for the CDW state since the collective 
modes are associated with poles of the density-density response function which can be derived in the time- 
dependent Hartrec-Fock approximation (TDHFA). Later, the TDHFA has been applied for anisotropic 
charge density wave (ACDW) states [28} and bubble states [29| to investigate the collective modes. 



3.1.2 vNL formalism for striped Hall, CDW, and bubble states 

Our vNL formalism for the CDW state is a similar one to the diagonalization technique developed by 
Yoshioka and Lee [261 ] . except for the use of the vNL basis instead of the eigenstate of the magnetic 
translation operator in one direction. Expanding the HF Hamiltonian by the vNL basis and adjusting 
the periodicity of the vNL to an integer multiple of the periodicity of the CDW, we can easily diagonalizc 
the HF Hamiltonian self-consistently. The essential difference is that a physical picture is more clear in 
the vNL formalism. In the vNL formalism, in contrast to other basis, momenta are clearly defined as 
eigenvalues of the commutative magnetic translation operators and the one-particle energy dispersion is 
obtained as a function of the momenta. The one-particle energy dispersion of the CDW states resembles 
that of noninteracting electrons in an external periodic potential [3(], [HI- I n the latter case, when the 
filling factor v (not the LL filling factor) is rational, i.e., for v = 4>q/BSq = q/p, where p and q have 
no factors in common and So = n with an electron density n is the unit-cell area of the crystal, the 
Landau level splits into p non-overlapping subbands. On the other hand, in the former case, while there 
is no periodic potential or external lattice structure, the von Neumann lattice plays a similar role and the 
one-particle energy dispersion of the CDW state has p bands for the partial LL filling factor v* = q/p. 
In both cases, the momenta arc defined in the Brillouin zone (or the magnetic Brillouin zone) and the 
one-particle energy dispersion has the p-band structure owing to the p-fold reduction of the Brillouin 
zone. In this sense, the physical picture is more clear for the vNL formalism. 

The vNL formalism was first developed for the striped Hall state by Ishikawa and others [32j, 1,0. 



Recently, this formalism has been developed for CDW states [12| and bubble states [35j. We review the 
details in the following sections. 
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Figure 3.1: Schematic view of the density profile of the striped Hall state. The density of the striped 
Hall state is uniform in one direction and periodic with the period ro in the other direction. If we take 
the vNL asymmetry parameter r s = ro, the HF Hamiltonian is diagonalizcd on the vNL basis. 

3.2 Striped Hall state 

3.2.1 Assumption of the unidirectional density 

Let us consider the case of the partial filling factor v* (0 < v* < 1) in the Ith LL. The striped Hall state 
is a unidirectional CDW state in the quantum Hall system, which has the following unidirectional density 
(Fig.EU): 

(pKx)} s tripe = E A K^)^(^ I 0)e i ( 2 ^ o )% (3.1) 

where ro is the period of the density in the x-direction, Ai(N x ) is an order parameter determined self- 
consistently, and Aj(0) = v* . Equation (|3.ip gives the following form of the mean value of the projected 
density operator, 

<P,(k)> stripe = Y, MN x )(2ir) 2 5(k x + ^-)5{k v ). (3.2) 

If we take the vNL asymmetry parameter r s = ro, Eq. (|3.2|) with k = (k x /r s , r s k y ) takes the simple form, 

(Mk))stri P e = ^AKAy^) 2 ^ + 2nN x )6(ky). (3.3) 

Nx 

3.2.2 Diagonalization, eigenstates and eigenvalues 

The HF Hamiltonian of the striped Hall state is already diagonalized on the vNL basis. Substitution of 
Eq. (|3.3p into the HF Hamiltonian Eq. (|2.41[) gives the HF Hamiltonian for the striped Hall state by 

^ ) F-st r i P e=/ ^^(pMp), (3.4) 
where e;(p) is a one-particle energy given by 

etip) = ef) + £ AKA>H^,0)(-l)^ e -^. (3.5) 
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v? F (o)/( q yh 



-1.25331 
-0.93999 
-0.80290 
-0.71968 



Table 3.1: Values of vf F (0) for each LL. 
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Figure 3.2: Fermi sea of the striped Hall state at v* . The occupied state is represented by the dark 
region. When the stripe faces the y-direction, the p^-direction of the Brillouin zone is fully occupied. In 
this case, the Fermi sea has the inter-LL energy gap in the p^-direction and is gapless in the p^-direction. 

In Eq. (|3.5p . ej ' is a uniform Fock energy given by ^*w ; HF (0). The values of ^ HF (0) are shown in Table 

ED 

For the ground state of the striped Hall state, the two-point function of the operator 6j(p) is given 
by [HS1 



<&i(p)6r(p'))stri P c = J>,J'0[cf - e i (p)](2^) 2 5 2 (p - p' - 2^N)e~^P< N \ (3.6) 

N 

where ep is a Fermi energy and 9 is a step function. 

3.2.3 Self-consistency condition and solution 

The self-consistent equation for Ai(N x ) is obtained by substitution of Eq. (|3.6p into the left hand side of 
Eq. (|3.3j) . Ai(N x ) = (~l) N:c sm{y* itN x ) / itN x is a solution of the self-consistent equation. This solution 
has the Fermi sea, \p y \ < nv* (shown in Fig. I3.2[) and gives the one-particle energy as 



(p) = 4 0) + E tn— ,o) sin( ^^ ^ 



ttN x 



(3.7) 



The HF energy per particle is given as a function of r a by 

C^Hg) stripe _ (^HF-stripc) stripe _ 1 (0)1 



F {1) (r \ 
-^stripeV 



2N, 



(0 



1 ,(0) , L V" *,.HF, 27r ^ / sm(v*wN x ) 



(3.8) 
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Figure 3.3: (a) 3D plot of the one-particle energy of the I = 2 striped Hall state at half-filling. When 
the stripe faces the y-direction, the one-particle energy is uniform in the p^-direction. (b) 2D plot of the 
one-particle energy as a function of p y . The Fermi velocity is logarithmically divergent. 
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I rf^ E stripc /(q 2 /l B ) 

1.636 -0.4331 

1 2.021 -0.3490 

2 2.474 -0.3074 

3 2.875 -0.2800 



Table 3.2: Minimum energy and corresponding parameter r s of the striped Hall states at v* = 1/2. 



P(X) 




Figure 3.4: Density profile of the I = 2 striped Hall state at half-filling. The uniform part po = v* is 
subtracted. The density of the striped Hall state is uniform in the y-dircction and periodic with a period 
r s in the x-direction. 

where is the total number of particles within the Zth LL. We determine the optimal value of r s by 
minimizing -Stripe ( r s) with respect to r s . The optimal value of r s and the minimum energy at each LL 
are shown in Table 13.21 (33| . 

For the optimal value of r s , the one-particle energy and the density profile of the striped Hall state at 
v* = 1/2 in the I = 2 LL are plotted in Fig. 13.31 and Fig. 13.41 respectively When the density is uniform 
in the y-direction, the one-particle energy is a function of only p y and uniform in the p^-dircction. The 
energy dispersion becomes gapless in the p y -dircction and has a inter-LL gap in the p^-dircction. The 
Fermi velocity vp = 9e;(p)/i9p; / |p B=TI y* is logarithmically divergent [H]. 



3.3 Charge density wave state 
3.3.1 Assumption of the CDW density 

Let us consider the case of v* in the Zth LL. Here, we consider only the rectangular CDW state for 
simplicity. The CDW state has the following periodic density (Fig. 13. 5[) , 

(p/(x))cdw = A/(QN)F ; (Q N )e-^ N ' x , (3.9) 

N 
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Figure 3.5: Schematic view of the density of the rectangular CDW state. The density is periodic with 
the periodicity r$ x and r$ y in the x-direction and y-direction, respectively. 

where Qn = (mQo x ,nQo y ) and N = (m,n) with integers m,n. Equation (|3.9p gives the following form 
of the mean value of the projected density operator, 

(pi(k))cDW = A;(QN)(2^) 2 <5 2 (k - Qn)- (3.10) 

N 

The reciprocal vector Rn' is given by Rn' = 2ir(m' /Q 0x , n'/Q 0y ) = (m' r 0x , n' r 0y ) , which satisfies 
Qn • Rn' = 2irM with an integer M . Since one electron in the CDW state with the density Eq. (|3.9[) 
occupies the area Sq = ro x ro y and one state in each LL occupies the area = a 2 with the vNL constant 
a = y h/eB, the partial filling factor v* is given by 

Number of electrons in the Zth LL a 2 

v = : — = . (3-11) 

Number of states in one LL roxroy 

Thus, the relation r$ x = a 2 /v*ro y holds. In particular, in the case of v* — N/M, this relation is given 
by tq x = M/Nro y , where we set a = 1 and the integers M,N have no factors in common. By using this 
relation, Qn is rewritten as 

„ , 27TTO TV . , 

Qn= ,2Trnr 0x — . 3.12 

ro x M 

Hence, (p/(k)) C DW is given by 

(p/(k)) C DW = E ^{ — ^nr^^fSikx - —)S(k y - 2 mx v*), (3.13) 

with v* = AI/N. If we take r s = r 0x (Fig. 13. 6[) . (p/(k))cDW is given by 
(pi(k)) C DW A;(^^, 2Trnr s v*)(2ir) 2 5(kx ~ 2Trm)5(k y - 2irnv*) 

oo M- l ON N 

= V V A l ( — : 2nr s (jM + k) — )(2TT) 2 S(kx-2nm)6(k y -2n(jM + k) — ), (3.14) 
^-^ r s M M 

m,j — — oo k— 



where n = jM + k with integers j, k (k = 0, 1, 2, . . . , M — 1). 
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electron 




Figure 3.6: The vNL unit cell and the CDW unit cell in the case of v* = 1/2. The thin lines represent 
the vNL. For v* = 1/2, the CDW unit cell is just twice as large as the vNL unit cell. If we take r s = ro x , 
the HF Hamiltonian can be block-diagonalizcd. 



3.3.2 Diagonalization, eigenstates and eigenvalues 

Substituting Eq. (|3.14|) into the HF Hamiltonian Eq. (|2.4ip . we can block-diagonalizcd the HF Hamilto- 
nian as follows. Under the assumption Eq. (|3.9p at v* = N/M, the HF Hamiltonian is written as 

(l) f d 2 k HF - 

a HF-CDW — J (2tt) 2 ' W 

x J2 E AK^^Trr^jM + A:) — )(27r) 2 J(fc x + 2™) ( 5(fc, + 2^(jM + A : ) — )p z (k), 

m,j — — oo k—0 

oo M-l AT 

=W>+ E E Ap F (m;,,^(-^,-2 7 rr s ( J A/ + fc)^), (3.15) 

m,j — — oo k—0 



where 

AHmjj.Aj^AK— ^^(jM + fc)^)^-— ,-2irr s (jM + k)^-), (3.16) 



except for m = j = k = and A} iF (0;0,0) = 0. Substituting Eq. (f2T40|l into Eq. ([3~T5]) and using the 
boundary condition of bi(p) given in Eq. (|2.29p . we rewrite the HF Hamiltonian as 



. , 2 oo M-l 

wSLodw - W> = / E E Af F (m;i,fc)e^W + W^ 

Jbz \ I m j—_ oa fc=0 

x bl(p)h(p x , Py + 2nk^)e im ^-^ N ^. (3.17) 
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Here, if we divide the interval of the p^-integral into M same intervals and reduce each interval to the 
range [0, 2tt/M] using the boundary condition Eq. (|2.29|) . H{|p_ CDW is given by 



n HF-CDW fc 2 Jv e — 



d Px r 2 * /M d P y 



2tt 



OO M-l 



2tt 



j7T(m+jN+m(j+{k-k')/M)N) 



vci ,j— — oo k,k'—0 



x b\(p x ,p y + 2^'^^^ + 2 7Tk^)e lm ^-^ N ^ + ^ mk ' N / M 



2- 



M-l 



2tt 



k,k'=0 



where 



N 

b(p,k) = b(p x ,p y + 2lT k—), 



OO OO 



A nF (p,k',k)= J2 E &HF(m;j,k-k>y^ m +3 N + 



m{j+{k+k')/M)N) Jmp y -i(jN)p<e 



(3.18) 

(3.19) 
(3.20) 



m— — oo j — — oo 



Hence, the HF Hamiltonian W^p_ CDW is block-diagonalized and the problem is reduced to the diago- 
nalization problem of the M x M Hcrmite matrix Ap F (p, fc', k). Note that the fact that (pi{r)) is real 
gives the relation A ; *(— Qjv) = A;(Qat), and by using this relation, it is proved that the M x M matrix 



a»*(p) with{AH P )} 



k'k 



Ahp(p, k',k) is Hcrmite. 



The Hermite matrix A; (p) is diagonalized using the unitary matrix U (p) as 



A(p) = C/t(p)AHp)[/(p) = 



/ co(p) 



V 



ei(p) 



(3.21) 



£M-l(p) / 



here, e s (p) represents the eigenvalue of the (s + l)th energy band, and the eigenvector for e s is given by 
{v s (p)} s / = U s ^ S '(p). RBZ represents the reduced Brillouin zone, < p x < 2ir, < p y < 2n/M. Hence, 



the HF Hamiltonian TL 



(I) 



HF-CDW 



is diagonalized as 



a HF-CDW fc i Jv e 



d 2 p 

ibz (2tt) 



-cJ(p)A(p)c ; (p) 



d 2 p 



M-l 



RBZ W ^ 



E e s (p)4(p) Cs (p), (3.22) 



where 



c;(p) = Mp),ci(p), • • • ,c M -i(p))* = C/ f (p)b ; (p) with 



M-l 



{b ; ( P )} fc = M P ,fc), Ca ( P ) = E ul k (pMp,k). 



(3.23) 



fe=0 



In the present case of v* — N/M, the lower N bands (from s = to s = N— 1) are fully occupied for the 
ground state. 



3.3.3 Self-consistency condition 

The HF Hamiltonian should be diagonalized self-consistently. The ground state is given by 



iV-l 



|n) = jv c l c t(p)|o), 



(3.24) 



s=0 pGRBZ 
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where N c is a normalization constant and |0) is a vacuum state in which the (I — l)th and lower Landau 
levels are fully occupied. Using this ground state, we obtain the self-consistency condition 

(n|p,(q)|fi) = J2 A i (QN)(27r) 2 ^ 2 (q - Q N ). (3.25) 

N 

In what follows, we derive the self-consistent equation for A(Qn) from Eq. (|3.25D . 

The field operators b](p) and bi(p) satisfy the anti-commutation relation Eq. (|2.30|) . For p and p' 
within the RBZ, Eq. ([230]) becomes 

{b\(p, k), b v {p', k')} = (2ir) 2 6 u ,5 k , k ,5 2 (p - p'). (3.26) 

For c|(p) and c s (p), the corresponding relation is given by 

{ct(p),<v(p')} = S s ,A^) 2 S 2 (p - p')- (3.27) 

Using Eq. ()2.40p . we rewrite the self-consistency condition Eq. (|3.25|) as 

M-i N 
(n\b}(p)h(p')\n) = £ F(p 1 k)Y,(^) 2 S(p x -p' x -2nN x )S(p y -p' y -2nN y -2nk—)e^^\ (3.28) 

fc=0 N 

where 

F(p,k)= Y Aj(— ,27rr a (jA/ + fe)— ) e -^{m+jN)-^m{j + k/M)N ^ (3 29) 

r s ' M 

m,j— — oo 

The left hand side of Eq. (|3~2"%1) is calculated using Eqs. (|3~2"3"1) and ([3~^7) . and the result is 

(Q|&}(p,fc)6*(p',0)|O) = r£ Ul^UoAp)^ (2n) 2 5 2 (p-p>), (3.30) 

for p and p' within the RBZ. From Eqs. (|3.25[) and (|3.30p . the self-consistent equation for Aj(Qjv) is 
given by 

A,(— ,27rr.(7M + = I ^2 ( J] ^ k (p)U , s (p)] jP*UN)-im Pve i« {m+j N)+i«m tf+fc/")". 

r s M J BZ (2tt) 2 \j^ Q s ' l J 

(3.31) 

If the right hand side of Eq. (|3.3ip coincides with the left hand side, the problem is self-consistcntly 
solved. 



3.3.4 Density profile and Energy dispersion of CDW states 

As in the case of the striped Hall state, the HF energy per particle is given as a function of r s by 

The optimal value of r s is determined by minimization of -EcdwO's) ^° Ts ' ^ nc optimal value of r s and 
the minimum energy per particle at several fillings in I = 0, 2 are shown in Table [3~"31 and Fig. 13.81 where 
cohesive energy, which is defined as (Total energy) — (Uniform Fock energy), is shown for the minimum 
energy. As for I = LL, the optimal value of r s is approximately given by y/l/i/*, which gives tq x = tq v . 
Hence, the CDW states in the I = LL become isotropic for the x and y-directions. On the other hand, 
as for the I = 2 LL, the optimal value of r s takes a different value from 1/v* especially near v* = 1/2, 
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which means that anisotropic states are energetically favorable at these fillings. Hence, the CDW states 
in the / = 2 LL become anisotropic near v* = 1/2. Since the original Hamiltonian has a rotational 
symmetry, the direction of the anisotropy can take any direction. In the present case of the rectangular 
vNL, two solutions whose anisotropy direction in density face the x or y-direction are degenerate. This 
degeneracy yields two optimal values of r s for anisotropic states. When we set the smaller one to 7g 
and the larger one to r s , they are related by 

r« = — Ly. (3.33) 

Only the smaller value ri^ is shown in Table 13.31 as the optimal value of r s for the CDW states in the 
I = 2 LL. 

The density profile of the CDW state is given by 

27rm JV 27rm TV 2nm N 

(Pi x) cdw = V Aj( ,2irr s n—)Fi{ , 2-KT s n— ) cos( x + 2Trr s n— y), 3.34 

z — ' r s M r s M r s M 

m,n— — oo 

where the x and y-inversion symmetry is assumed. In Figs. 13.91 and I3.10[ the density profile is plotted 
at several fillings in / = 0, 2 LLs for the optimal value of r s . As for the I = CDW density, the density 
becomes isotropic for all partial fillings. On the other hand, as for the I = 2 CDW density, the density 
becomes highly anisotropic near v* = 1/2, and as the partial filling decreases, the density approaches the 
isotropic one gradually. This difference between the lowest LL and the higher LLs is mainly due to the 
difference of the form of the HF potential w ; HF (k). We refer to the direction with a narrower periodicity 
in density as "stripe direction" in this thesis. 

In the present vNL formalism, the energy dispersion of the CDW state is obtained as a function 
of momenta p x and p y which are conjugate to the vNL. As for the I = LL, the energy dispersion is 
isotropic in the p x and p y -directions and the band width becomes narrower as the density becomes lower. 
In Fig. I3.11[ the energy dispersion at p y = is plotted as a function of p x . As for the I = 2 LL, the 
energy dispersion becomes almost flat in the direction perpendicular to the stripe direction. In Fig. I3.12[ 
the energy dispersion at p y = is plotted as a function of p x . In both cases, there is a large energy gap 
between the filled bands and the empty bands. This large gap coincides with the previously obtained 
results by Yosioka and Lee [26[ . 



3.4 Bubble state 

3.4.1 Assumption of the bubble density 

The bubble state is a CDW state with two or more electrons per unit cell. We call the bubble state 
with TV electrons per unit cell the "TV -electron bubble state" . Let us consider the TV-electron rectangular 
bubble state at v* = l/M in the Ith LL, with a period ro x and ro y in the x and y direction, respectively. 
In this case, since TV electrons are in the area of the unit cell Sq = ro x ro y and the density of states is 
given by 1/a 2 , the partial filling factor is given by 

^x-^y 7\7" q 

„ Number of electrons in the Zth LL r xi- 0y a ^ 

Number of states in one LL L x L y ro x ro y ' 

where L x and L y are the size of the system in the x and y directions, respectively. In particular, in the 
case of v* = l/M, this relation becomes r$ x = NM/ro y , where we set a = 1. By using this relation, the 
mean value of the projected density operator of the TV-electron bubble state is given by 

°° 27rm 1 1 

(pz(k))bubbic= Yl M^,2™r s — )(2ir) 2 5(k x -2Trm)5(k y -2irn—) 7 (3.36) 
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where we set r s = To- Hence, the diagonalization procedure of the HF Hamiltonian is the same as in the 
case of the CDW state at v* = 1/K with K = NM. 



3.4.2 Self-consistency condition and solution 

The important difference from the CDW state is a self-consistency condition. The one-particle energy of 
the A'-electron bubble state at v* = 1/M has K = NM energy bands. Since the filling factor is 1/M, 
the lowest N bands are fully occupied for the ground state. As a result, the self-consistency condition is 
given by 

A,( — ,2^(n + ^)r s ) = / — V U* s (p)C/ , s (p) (-l)™+»+™« e ^»-^— -™07*). (3.37) 
r s & J BZ (2tt) \ s=Q ■ J 

For the I = LL, no self-consistent bubble solution has been obtained. For the / = 2 LL, a self- 
consistent bubble solution is obtained for v* = 1/4. The obtained solution is a 2-elcctron bubble state. 
The optimal value of r s is r s = 2.8 and the corresponding cohesive energy is —0.1576 (q 2 /Ib) which is 
slightly lower than that of the CDW state at the same filling. Hence, for the v* = 1/4 in the / = 2 
LL, the 2-clectron bubble state is more stable. The density profile of this state is shown in Fig. 13.131 
which is almost isotropic. The energy dispersion consists of eight bands and the lower two bands arc fully 
occupied. As in the case of the CDW state, there is a large gap between the filled bands and the empty 
bands. 

Since electrons of this 2-electron bubble state are well-localized in density, it is expected that the 
state is easily pinned by impurities, which results in a quantization of the Hall conductivity owing to its 
collective insulating property. Indeed, this quantization has been observed around v* = 1/4,3/4 in the 
/ = 2 LL [13, EH, and ^ 

is believed that this observed quantization is an evidence of the bubble state. 



3.4.3 Relation between the iV-electron bubble state and the CDW state at 

v* = 1/M 

The diagonalization formalism for the A^-electron bubble state at v* = 1 /M includes the CDW solution 
at v* = 1/M. To make this point clear, let us compare the mean value of the density of the CDW state 
and the bubble state. For the CDW state at v* = 1/M, the mean values of the density operators are 
given by 

<Mx))cdw= E AP DW (^,2.r c n-i)F i (^,2^ c ni-) e -( 2 -/-^-( 2 ^/^^, (3.38) 



ra,n— — 00 
00 



<Mk)>CDW= £ A?™(^,?^)(2nfS(k x -2 7 rrn)S(k y ~2nnj-), (3.39) 
^-^ rc M M 



m 1 n——co 



and for the A^-electron bubble state at v* = 1/M, they are given by 

<Mx))bubb lc = £ A™ e (^,27rr B nl)FK^,2 7 rr B ni)e- i ( 2 ^B)x-i(2 7rrBra /K), ) (34Q) 



m,n— — 00 
00 



(Pl(k))bubble = E ^ uhhl %^,^^)(^) 2 S(k x -27rm)S(k y -27rn^), (3.41) 



m,n— — 00 



with K = MN, where r s = rc for the CDW state and r s = r-g for the bubble state. If we set re = 
Nr c , A\ nhhlc {2TTNm/r B ,2irr B n/K) = Af DW (27rm/r c , 2irr c n/M ) and other Af ubblc s is equal to zero, 
(pi(x))bubbic becomes the same as (p;(x))cdw- Hence, the HF Hamiltonian with the assumption (|3.36[) 
becomes the same as the HF Hamiltonian with the assumption (|3.14|) . 
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Figure 3.7: Schematic view of the density of the 2-electron bubble state. 
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Table 3.3: Minimum energy and corresponding parameter r s of the CDW states at several fillings in the 
I = 0, 2 LL. As for the minimum energy, the uniform Fock energy is subtracted. The optimal value of r s 
for the / = 2 LL takes two different values which give the same energy. Only the smaller one is shown in 
this table. These values are plotted in Fig. 13.81 




Figure 3.8: The optimal value of r s is plotted for Z = 2 CDW states. Cohesive energy is plotted for I = 
CDW states with the optimal value of r a . 
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Figure 3.9: 2D plot of the density of the I = CDW state. As the filling decreases, the distance between 
localized electrons increases isotropically. 
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Figure 3.10: 2D plot of the density of the I = 2 CDW state. Near v* = 1/2, the density is highly 
anisotropic. As v* decreases, the density becomes isotropic gradually. 
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Figure 3.11: 2D plot of the one-particle energy dispersion of the / = CDW state. The energy dispersion 
at p y = is plotted as a function of p x . There is a large gap between the filled bands and the empty 
bands. 
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Figure 3.12: 2D plot of the one-particle energy dispersion of the I = 2 CDW state. The energy dispersion 
at p y = is plotted as a function of p x . There is a large gap between the filled bands and the empty 
bands. 
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Figure 3.13: 2D plot of the density of the I = 2 2-electron bubble state at v* = 1/4. The density is 
almost isotropic in the x and y directions. 
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Figure 3.14: 2D plot of the one-particle energy dispersion of the I = 2 2-electron bubble state at v* = 1/4. 
The energy dispersion at p y = is plotted as a function of p x . There is a large gap between the filled 
bands and the empty bands. 
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Chapter 4 

Effects of an injected current on 
highly anisotropic states 



In this chapter, we study effects of an injected current on highly anisotropic states which have been 
discovered around half-filled third and higher Landau levels (LLs) in experiments with ultra-high mobility 
samples at low temperatures. While several states have been proposed to theoretically explain the 
experimental results so far, it is still an open problem which state is realized in the experiments. Most 
previous studies have been done within the Hartree-Fock (HF) approximation with no injected current. 
In experiments of anisotropic states, however, a current is injected to investigate properties of the system. 
This effect has not been taken into account in the previous studies. In this chapter, we focus on two 
different anisotropic HF states, i.e., a striped Hall state and an anisotropic charge density wave (ACDW) 
state, and investigate the effect of the injected current on these two HF states using response functions 
against the injected current. The calculations are done by means of a von Neumann lattice (vNL) 
formalism developed in Chapter [2] With no injected current, the ACDW state has a lower energy. We 
find that the striped Hall state becomes lower energy state when the injected current exceeds a critical 
value. The critical value is estimated at about 0.04-0.05 nA, which is much smaller than the current used 
in experiments of anisotropic states. 



4.1 Introduction 



In the quantum Hall system, half-filled states at each Landau level (LL) exhibit much attractive features. 
Around the half- filled lowest LL, isotropic compressible states have been observed [3(1 which are 
widely believed to be the Fermi liquid of composite fermions. Around the half- filled second LL, the 5/2 
fractionally quantized Hall conductance has been observed [3^|. The p-wave Cooper pairing state of 
composite fermions, which is called the Pfaffian state, has been proposed to explain this state [39], |40| . 
Around half-filled third and higher LLs, highly anisotropic states, which have extremely anisotropic 
longitudinal resistivities and un-quantized Hall resistivities, have been found in ultra-high mobility sam- 
ples at low temperatures M , 1 1 011 - M uch theoretical work has been done to study the anisotropic states 
H, [H, EH, SI, 0, [44l.[45Tl46l.l47l. 0, [isl and several states have been proposed so far. In this chapter, 
we focus on two different HF states among them, i.e., a striped Hall state [1, [4l[ and an anisotropic 
charge density wave (ACDW) state [28| . 

The striped Hall state is a unidirectional charge density wave state which is a gapless state with an 
anisotropic Fermi surface [32I [33[ . The Fermi surface has an energy gap in one direction and is gapless 
in the other direction (Figs. l3~2l and l4~Tj) . The ACDW state is a CDW state which has a similar periodic 
density in a direction perpendicular to stripes, and in addition, has a density modulation along stripes 
(Fig. 14. 2[ ). The density modulation along stripes results in energy gaps in both directions. These features 
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suggest that the anisotropic state found in experiments is the striped Hall state since the anisotropic 
longitudinal resistivity and the un-quantized Hall resistivity are naturally explained by the anisotropic 
Fermi surface [3^ , |33| , while it is difficult to explain these experimental features with the ACDW state 
because of the energy gap. However, it has been pointed out that the striped Hall state is unstable within 
the HF approximation to formation of modulations along stripes so that the ACDW state is the lower 
energy state (l^. This has been an enigma as to the anisotropic states. 




Figure 4.1: Left: Density modulation of the striped Hall state at the half-filled third LL. The uniform part 
is subtracted. The density is periodic in the cc-direction and uniform in the y-direction (stripe direction). 
The injected current flows easily in the stripe direction. Right: Energy spectrum of the striped Hall state 
at the same filling. It is uniform in the p^-direction. The BZ is \pi\ < n and Ep is a Fermi energy. 
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Figure 4.2: Left: Density modulation of the ACDW state at the half-filled third LL. The uniform part 
is subtracted. The density is modulated along stripes. Right: Energy spectrum of the ACDW state at 
the same filling. It depends on p x weakly and the spectrum at p x = is plotted. The BZ is reduced to 
\p y \ < ir/2 owing to the modulations along stripes and the two bands are formed. 

In the experiments of the anisotropic states, current is injected. This effect has not been taken into 
account in the previous calculations of the total energy. About two decades ago, MacDonald et al. have 
studied the effect of the injected current on the integer quantum Hall system [50(. They calculated the 
current and charge distributions and found that charges accumulate around both edges of the sample 
with the opposite sign, as expected from the classical Hall effect 5^, 5l|, |B2, H3] • The charge accumulation 



causes the energy enhancement via the Coulomb interaction between charged particles. The same type 
of energy corrections may exist even in the present highly correlated quantum Hall states. However, the 
effect of the injected current on the anisotropic state has not been studied. 

In this chapter, we calculate the correlation energies of the striped Hall state and the ACDW state 
in the system with the injected current flowing in the direction along stripes, and with no impurities and 



4.2. EXPERIMENTS OF HIGHLY ANISOTROPIC STATES 



37 




Figure 4.3: Longitudinal resistivity p xx in sample A at temperature T = 25 mK for two perpendicular 
directions of the current flow through the sample. The two traces result from simply change of the 
direction of current through the sample; the sample itself is not rotated. The diagrams in the figure 
depict the difference between two configurations. From Ref. 0. 



no metallic contacts. It is important to know if the ACDW state has a lower energy even in the system 
with the injected current. The density modulation along stripes shown in Fig. 14.21 suggest that the effect 
of the injected current becomes larger in the ACDW state. To verify whether this naive expectation 
is true or not, the dependence of the correlation energies on the injected current is studied in detail. 
Effects of impurities and metallic contacts are ignored in our calculations of correlation energies since 
these effects are expected to be small in experiments of anisotropic states in ultra-high mobility samples 
and are outside the scope of this work. Effects of the injected current are investigated using response 
functions for electromagnetic fields. The current and charge distributions are determined and energies of 
the two states are calculated from these distributions. It is found that the energy of the ACDW state 
increases faster than that of the striped Hall state as the injected current increases. Hence, the striped 
Hall state becomes the lower energy state when the current exceeds the critical value. The critical value is 
estimated at about 0.04 - 0.05 nA, which is much smaller than the current used in the experiments. Our 
result suggests that the anisotropic states observed in experiments are the striped Hall states. Hence, 
the enigma as to the anisotropic states is resolved. 

This chapter is organized as follows. In Sec. 14.21 we review the experiments of the anisotropic states 
and discuss the problem briefly. In Sec. 14. 3[ properties of the striped Hall state and the ACDW state are 
discussed within the HF approximation with no injected current. In Sec. 14. 4[ electromagnetic response 
functions of the two HF states are calculated in the long wavelength limit. Using these response functions, 
we determine the current and charge distributions and calculate the energy corrections due to injected 
currents in Sec. 14.51 A summary is given in Sec. 14.61 



4.2 Experiments of highly anisotropic states 

In 1999, working with 2D ultra-high mobility samples at low temperatures, Lilly et al. @ and Du 
et al. [101 ] have reported strong anisotropy in longitudinal resistivities p xx around v* = 1/2 in I > 2 
LLs. The experimental result by Lilly et al. is shown in Fig. 14.31 In this figure, the sample A is a 
GaAs/AlGaAs hcterojunction grown by molecular beam epitaxy, where the density of the sample is close 
to n s = 2.67 x 10 11 cm -2 , the low temperature mobility is \i > 9 x 10 6 cm 2 /IAs, and the sample size is 
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a 5 x 5 mm square. The electrical transport measurements have been performed using 2-20 nA, 13 Hz 
excitation. For more details about experimental parameters used by Lilly et al. , see Ref. We refer 
to these experimental parameters as Lilly's parameters in this thesis. 

In this experiment, longitudinal resistivities p xx and Hall resistivities p xy have been measured for two 
different configurations. For the configuration shown in the left inset in Fig. 14.31 p xx is almost zero around 
half-fillings in / > 2 LLs, such as v = 9/2, 11/2, 13/2. On the other hand, for the other configuration 
shown in the right inset in Fig. l4.3( p xx revels very large values around the same fillings, where the sample 
itself is not rotated and only the direction of current is changed. Here, no quantized plateau in the Hall 
resistivity p xy have been observed around these fillings. These experimental results indicate that the 
ground state around these fillings is not an isotropic state such as a fractional quantum Hall state but an 
anisotropic state. Since the original Hamiltonian of the quantum Hall system has a rotational symmetry, 
it is considered that the anisotropic state is energetically favorable and a spontaneous symmetry breaking 
occurs with respect to the direction of the anisotropy which would be determined by some small residual 
inhomogeneity of the sample. 

4.3 Highly anisotropic states in the Hartree-Fock theory 

In this section, the striped Hall state and the ACDW state at the partial filling factor v* = 1/2 in higher 
LLs are investigated. Calculations are performed by means of a vNL basis within the HF approximation 
with no injected current. As has been mentioned in Chapter [21 the vNL basis is a suitable basis to study 
spatially periodic states. 



The striped Hall state is obtained as a HF solution under the assumption of the unidirectional charge 
density wave Eq. (|3.1|) . The solutions at v* = 1/2 are given in subsection 13.2.31 in Chapter [3] in detail. 
Let us consider the striped Hall state in which stripes face the y-dircction. In this case, the anisotropic 
Fermi surface of the striped Hall state has the inter-LL energy gap in the p^-direction and is gaplcss in 
the Pj,-direction. When a current flows in the y-direction, an electric field is generated in the x-direction 
owing to the Hall effect so that the Fermi sea slides in the p^-direction. Since the Brillouin zone in the 
p^-direction is fully occupied, there is no dissipation in this case and the longitudinal resistivity in the 
y-direction is expected to be zero. On the other hand, when the current flows in the x-direction, the 
Fermi sea slides in the py-direction so that the energy dissipation occurs owing to the gapless structure 
of the Fermi surface in this direction, which would result in a large longitudinal resistivity. Hence, it is 
expected that the striped Hall state revels the highly anisotropic longitudinal resistivity, which coincides 
with the experimental results. 

4.3.2 Anisotropic charge density wave state at v* = 1/2 

CDW states studied in Sec. 13.31 in Chapter [3] become highly anisotropic around v* = 1/2 in higher LLs. 
We call these highly anisotropic states at v* = 1/2 "anisotropic charge density wave (ACDW) states" 
and the direction with a narrower periodicity in density "stripe direction" in this thesis. In the vNL 
formalism, the HF Hamiltonian of the ACDW state can be rewritten in the form of a 2 x 2 matrix. 
Therefore, eigenvalues and eigenstates for the ACDW state can be easily obtained in a more analytical 
form as a function of order parameters. In what follows, we derive them and discuss the features of the 
self-consistent ACDW solution. 

Let us start by the HF Hamiltonian of the CDW state at v* = 1/2. From Eq. (|3.15[) , this is given by 



4.3.1 Striped Hall state at v* = 1/2 




(4.1) 
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Figure 4.4: One-particle energy of the / = 2 ACDW state at half-filling. The uniform Fock energy is 
subtracted. Two bands are formed and the lower band is completely filled. 

h,(v)-( bl{ P x ^ \ D,(u) = ( A(p) B(p) ^ (4 2) 

where e, is a uniform Fock energy given by ^*w ; HF (0), the momentum integration is performed over the 
reduced Brillouin zone (RBZ), \p x \ < tt and \p y \ < 7r/2, Di(p) is a 2 x 2 Hermite matrix, and ^4(p), Bp 
are given by 

A(p) = ]T vJ lF ('—^-, 2wNyr 3 )Ai(———, 2irN y r s )e iv ( Nx+Ny+NxNy ^~ ipxNv+iPvNx 

2ttN x _ /flAr , -i \ \ a ^TriVa _, nAr , ^ , Jir(N m +N s +N m (Ny+l/2))-ip x Ny+ip y N x 



r 

N 



B{p) =J2vr(^M2N v + l)r.)M^MW y + l)r 

(4.3) 



The Hamiltonian Eq. (|4.1[1 can be diagonalizcd at each momentum just by unitary transformation of the 
field operator b;(p). In the present case of v* = 1/2, the Brillouin zone is reduced to the half size of 
the original domain and two energy bands are formed (Fig. 14. 4[) . Di (p) is diagonalized using the unitary 
matrix U(p) as 

UHp)D l (p)U(p) = ( £ ^ P) £ ° {p) ) . (4.4) 

where e+(p) and e_(p) represent the upper energy band and the lower energy band respectively. If we 
assume the x- and y-inversion symmetries for the density, the order parameters become real and have 
the property Ai(Q x ,Q y ) = Ai(Q x ,—Q y ) — Ai(—Q x ,Q y ). The self-consistent solution is available only 
when the two energy bands are symmetric with respect to the energy gap, i.e., e+(p) = — e_(p) = e(p), 
as expected from the particle-hole symmetry of the original Hamiltonian. This gives Tr2 X 2-Di(p) = 0? 
where T~2 X 2 denotes the trace with respect to the 2x2 matrix indices, and A(p x ,p y + tt) = —A(p). In 
this case, U (p) and e(p) are given by 

/ B(p) B(p) \ 

U(P) = t+ N p )-%p) ^Hp) . <P) = V(MpW + \B( P )\2, (4.5) 

\ JV+(p) JV_(p) / 

where N±(p) — 2e±(p)(e±(p) — ^-(p))- Using the basis c;(p) = ?/T(p)b;(p), we obtain the HF Hamilto- 
nian of the ACDW state 

<-«cow + / ,0c»(p, ( e _» p) ) c«(p), (4.6) 
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I rf CDW E ACBW /(g 2 /l B ) A Gap /(g 2 /Z B ) 

V2 -0.4436 0.3292 

1 1.02, 1.96 -0.3583 0.3077 

2 0.82, 2.44 -0.3097 0.2470 

3 0.70, 2.86 -0.2814 0.1967 



Table 4.1: Minimum energy and corresponding parameter r s of the ACDW states at v* = 1/2. Ac a p is 
a magnitude of the energy gap. 



where 

and A;(Qjv) is determined as a numerical solution of the self-consistent equation Eq. (|3.31[) . 

The HF energy of the ACDW state depends on the asymmetry parameter r s . The optimal value of r s , 
the HF energy per particle, and the magnitude of the energy gap are given in Table [4TT1 in which there 
are two values of r s at each LL owing to the 7r/2-rotational symmetry. The magnitude of the energy gap 
is estimated to be of the order of 10 K for a few tesla. Therefore, the extremely anisotropic longitudinal 
resistivities are not expected at tens of milli Kelvin and it is difficult to explain the experiments with the 
ACDW state. On the other hand, the HF energy of the ACDW state is slightly lower than that of the 
striped Hall state at each LL. This was one of the remaining issues for the anisotropic states. 



4.4 Response functions 

In this section, electromagnetic response functions of the two HF states are calculated in the long 
wavelength limit. We consider the quantum Hall system with an infinitesimal external gauge field 
a^x) = (ao(x),—a(x)) and calculate the response functions of the striped Hall state and the ACDW 
state. These response functions will be used to determine current and charge distributions and energy 
corrections due to an injected current in Sec. 14.51 



4.4.1 Response function of the striped Hall state 

The Hamiltonian in the quantum Hall system with a M (x) is given by 

H=J d 2 x*t(x) ( HV + eA 2 ( ^ } + ga(x))2 - eao(x)) *(x) + i J d 2 xd 2 x' : />(x)V(x - x')p(x') :, 

(4.8) 

where V(x) = q 2 /|x|. We project the Coulomb interaction part to each LL and apply the HF approxi- 
mation to the projected Coulomb interaction. Then, by using the vNL basis, the Hamiltonian in the HF 
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approximation is given by 



bz (2^) 2 

/ (0E<,(kK(k) / M ^Cp )6 ,(p-k)e-(V^-W 
/ 7^ E • + k') |^I(pMp _ k _ ^ )e -(i/^ + K)^ v -K-K) 

E / 702< F (k)(Pi(-k))w(k), (4.9) 



where i2 (k) is defined by (see Appendix IA.4I) 

//U(k) H/d^V^}^) (4.10) 

in which i/ 1 = (1, — lu c t], w c £) is the electron velocity. Repeated Greek indices \i and v are summed for 0, 
1, 2 in throughout the following calculations. The action is given by 

S[a li ,b,b*] = j dt[J^^bj(p,t)(id t + p F )b l (p,t)-H(t)\, (4.11) 

where H{t) is the Heisenberg representation of H . 

Let us concentrate on the striped Hall state with the stripe direction facing the y-direction at half- 
filling. Substituting Eq. (|3.3p into Eq. ()4.11j) . we obtain the action of the striped Hall state given by 

SHp[o M> 6,6t] = W ^it (p) L _ (: l{p)) 6 u ,(2nf6 3 (p- P >) - U^(p,p') - E/^'W)1 htf), 

(4.12) 

where 

ttf/W) =-J2 e //>(P P' 2 7 rN)/ l (p + p',NK(po -K,p - P ' - 2TrN)e-W 4 "fo'-M<n> + ' / *\ 

N 

d 3 fc 



^ i ' ) (p>p / )=E / (^^'(p-p'-^Mp + p'-N) 



N 

x a(fc , k) • a( Po - p' - k , p - p' - k - 27rN)e- (4/47r)(p ^^ )(p " +p » ) 

/i(p,N) = f y -i) N ^+ N v+ N ^ N v e -( i / 2 )p^ N v+( i / 2 )Pv N ^ . (4-13) 

Here, p denotes (po, p), £z(p) = Ei + e;(p) — £f, and (p,p') and [/a 2 '' \p,p') are the first order term 

and the second order term with respect to a M (p), respectively. 

The partition function Z[a] is calculated using path integrals by 

Z[a„] = j vtfVbe iSm ^> h ^ ] 

= J X>6tx) foe -(-»)Ef' t (s- 1 l-C/ ai -C/ a2 )b = e T l -log[(- t )g- 1 ] e Tr\og[l-gU ai -gU a2 ]^ 

where the power of the exponent is expressed in the matrix representation in momentum space and Tr 
denotes the trace of the momentum indices and the LL indices. gi{p) is the Green function given by 

,* 0(6 (p)) 0(-6(p)) r ,,~ 

Po-&(p)+™ Po-Cz(p)-* 5 
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where S is an infinitesimal positive constant. The effective action S c ff[^] is defined as S^ffl^i] = 
— ilogZ[a AI ], which consists of the non-perturbed part Sq = Tr log[(— and the correction part 



(4.16) 



due to the external gauge field AS c s[afj]. AS c g[afj] is given by (Fig. 

AScff [a„] = A5i[a^] + AS 2 [a t ,} + AS 3 [ aA1 ] + 0[a% 
where 

ASAa,] =VTt\gU ai ] = / S^giW^fap), 

x JBZ \ Z7T ) 

AS 2 K] =iTr[gU a2 ] = iJ2 I S^9i(p)U^ l Hp,p), 

AS,[a,] =lTr[gU ai gU ai ] = \ ^^g^u^p^g.^U^ 1 ^' ,p). (4-17) 

Substituting the expressions for <#(p), U^{ 1 (p,p') and f/^*' (p,p') into Eq. (|4.17p . we obtain 




Figure 4.5: (a), (b) and (c) are Feynman diagrams for ASx, AS2 and AS3, respectively. 



ASWtK] =(eio)ao(0, 0) + ]T e/^, (-27rJV x /r„ 0) 



AT* 



sin(p F iV a ) ci7rjVa 

7rAL 



0^(0,-2^/^,0) 



5 / ^^«M(Po,p)^ /1I/ (p,N)a,(-p 0j -p-27rN), 



(4.18) 



where Zq represents the uppermost partially filled LL and N = (N x /r s ,r s N y ). K flv (p,'N) is a response 
function given by 



K»»(p, N) = <?f?,v (P)fv,i (-P - 2^N)/ M , (p , p, N)fc(p, N) 



/,/' 



2tt 



N 



*oOn,o + /i ,i (-27rJN)(-l) dAT yi0 



(^,1^1+^,2^,2), (4.19) 



and Ij^' (po, P, N) is the loop integral in Fig. 14.51 (c), which is given by 



i/,j'(po,p,N) 



BZ 



(2tt)' 



K*(p + p')M6'(p')) 



Po + 6'(P') - 6(P + P') + W Po + fr(P') - 6(P + P') - iS 



(4.20) 
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in which the p' integral has been performed. In Eq. (|4.19j) . the first term and the second term come from 
AS^a^] and AS^a^], respectively, and the second term is canceled with the p = part of the first term, 
as expected from gauge invariance. Hence, K l ^(p = 0, N) = for i = 1 , 2. 

In the long wavelength limit, the largest contribution in the response function comes from the N = 
part. In the case of po = p y = and p x — > 0, which is used in the next section, the largest contribution 
comes from the lowest order term in Kq (p x ) = K^ v '(0,p x , 0) with respect to p x . Expanding Kq \p x ) up 
to the lowest order, we obtain the response functions in the long wavelength limit given by 

K (P*) = ^Pt, K°y( Px ) = -ia^p x , 

Kf( Px ) =ia%) Px , Kf{p x ) = a K ^ cV l, (4.21) 

where o xy — e 2 v/2-K and ax = e 2 w c identified as the Hall conductance 

since if we consider a static homogeneous electric field in the x-direction generated by the gauge field 
Oq x (x) = xE x , then the electric current in the y-direction (j y (x)) is given in the long wavelength limit 
by (j y (x)) = 5AS c g/5a y (x) = Kq° (d x )ao(x) = —a xv 'E x , where the response function transformed in the 
coordinate space is used. The longitudinal resistivity becomes zero in the present calculation since the 
impurity potential is not included. If impurities are added, it is expected that the longitudinal resistivity 
becomes zero in one direction and finite in the other direction owing to the anisotropic Fermi surface. 



4.4.2 Response function of the anisotropic charge density wave state 

The action of the ACDW state at half-filling is given by 



where V^' 1 ' [p;p') with j = 1, 2 is a 2 x 2 matrix given by 



E 



d 3 P d 3 p' t 

RBZ { Z7r > 



G- 1 ( P )S u ,(27TfSHp~p / ) - V«/\p;p') - V^(p; P ')] Cl ,(p'), 

(4.22) 



K ( f W) =i^p)uifWMp), 



u^'W) 



Ui 1 / \p;Po,p x ,p' y + 7r) 



Ua/ \po,p Xl Py + tt;p') Ui 1 / \po,p x ,p y + TT;p' Ql p' x ,p' y + tt) 



(4.23) 



and Gi (p) is a 2 x 2 matrix Green function given by 

Gi( P ) = 



9iip) 
grip) 

gr(pn 

gtip) 1 



for I = Iq 

for I < lo 
for I > Z 



(4.24) 



In Eq. (|4.24p . g^{p) = l/(po — £,t(p) ^ i s a one-particle Green function for the upper or lower band, 
1 is a 2 x 2 unit matrix, and ^(p) = Ei a + e£°^(p) - e F . 

The partition function Z[a^\ is calculated using path integrals by 



p c t : p ce i'S , HF[<V,c,ct] = e T>log[(-i)Gt- 1 ] e 'Erlog[l-GVa 1 -G!V 03 ] < 



(4.25) 



The correction part of the effective action AS c ff [a M ] is given by 



AS eff [a M ] = A5i[a M ] + AS 2 [a M ] + AS 3 [a^} + 0[a% 



(4.26) 
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where 



AS 1 ! [a M ] =iTr[GV ai ] = iY i f 

i ^ 

A^K] =iTr[GV a2 ]=iJ2 



d 3 p 

;Z W 
d 3 p 



RBZ 



(2tt) 



A5 3 [a M ] =>[GV ai GV ai ] = ^ 



/./' 



RBZ 



■&ax 2 [G? I (p)V a ?> ,) (p;p)] ) 

Tr 2x2 [G ; (p)K ( 2 M) (p;p)], 

[G l (p)VW')(p;i/)Gr Jf (p / )Vi 1 , '' , )(p';p)]. (4.27) 



d 3 pd 3 p' 



(27T) e 



Tr 2 



Here, Tr denotes the trace with respect to the momentum indices, the LL indices and the 2x2 matrix 
indices, and Tr 2x2 denotes the trace with respect to only the 2x2 matrix indices. Substituting the 
expressions for Gi{p), Va[' 1 \p;p'), Va/ \p;p') into Eq. (|4.27|) . we obtain the following expression for the 
N = part of ASefffa/J as m the case of the striped Hall state, 



AS'cffK] =e(Z o + ^)ao(0,O) 



d 3 p 

(27T) 



3 a ll (p ,p)K^ l '(p)a l/ (-pQ, -p), 



(4.28) 



where Kq V (p) is given by 



EE 



;<i £'>(o 
1 



uS> - 1) 



+ 



E+E 

;<z i>io. 
d 2 p' 



{/^(p)/^(-p) + /^(p)/m'(-p))} 



{/CJp)/^(-p) + /£>)/mo(-p)} 



w c (/ - 
1 



rbz (27r)2p -( e (p') + e( P + p')) 

A(p + p')^(p') + Rc(S(p + p')B* (p^e-^/ 2 ^* 



e(p + p')e(p') 



1 e 2 

Z Z7T 



(4.29) 



up to 0(e(p)/ujc)- In Eq. (|4.29p . the last term is canceled with the p = term of the second term as 
expected from gauge invariance again. Hence. A7,''(0) = for i = 1.2. 

In the long wavelength limit po = p y = and p x — > 0, the largest contribution in the response function 
comes from the lowest order term with respect to p x . The expressions of K v (p x ) and Kq (p x ) become 
the same as in the case of the striped Hall state. The expression of Kq°(p x ) becomes slightly different by 
the correction from the intra-LL effect at the uppermost partially-filled LL. For the striped Hall state, the 
one-particle energy shown in Fig. 13.31 has the inter-LL energy gap in the p^-direction, and the inter-LL 
effect gives the response functions given in Eq. (|4.19p . For the ACDW state, the one-particle energy 
shown in Fig. 14.41 has the intra-LL energy gap in the p^-direction as well as the inter-LL energy gap. 
While the inter-LL effect gives the same expression of the response function as that of the striped Hall 
state, the intra-LL effect causes some corrections to the response function. Including these corrections, 
we obtain K®°(p x ) of the ACDW state in the long wavelength limit given by (see Appendix IB. 2|) 



K»»( Px ) = -(l + 



V UJ r . 



(4.30) 



where a = y/2ii/eB is used explicitly in order to compare the theoretical results with experimental data. 
The value of (3 at each LL is shown in Table FOl Note that the unit of (3 is (tesla) -1 / 2 . 
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I parallel perpendicular 

0.497 0.497 

1 0.386 0.824 

2 0.472 1.044 

3 0.581 1.171 



Table 4.2: The value of (3 at each LL. parallel is the case in which the stripe direction faces the y-direction. 
perpendicular is the case in which the stripe direction faces the ^-direction. The unit of /3 is (tcsla) -1 / 2 . 

The Hall conductance is given by ai u J = e 2 i//27r, as in the case of the striped Hall state. The 
longitudinal resistivity becomes zero in the present calculation since the impurity potential is not included. 
However it is expected that the longitudinal resistivity remains zero even in the system with impurities 
because of the energy gaps. 



4.5 Energy corrections due to injected currents 

In this section, we consider the quantum Hall system with an injected electric current and investigate the 
current effect on the striped Hall state and the ACDW state. Effects of impurities and metallic contacts 
are ignored in our calculations. For the striped Hall state, we only consider the current parallel to the 
stripe direction since in this case, the current effect can be estimated with no ambiguity even in the 
system with impurities. When the current flows in the stripe direction, charges accumulate around both 
edges of the sample in the direction perpendicular to the stripe direction (perpendicular direction), as we 
will see later, and the electric field generates in the perpendicular direction. In this case, the impurity 
effect is negligible since the Fermi surface has the inter-LL energy gap in the perpendicular direction. 
On the other hand, when the current flows in the perpendicular direction, the impurity effect becomes 
relevant since the electric field generates in the stripe direction while the Fermi surface is gapless in this 
direction. The current effect in this case is nontrivial and will be studied in future work. 

In the system with an injected current, it is naively expected that the current flow causes the plus 
and minus charge accumulation at both edges of the sample with the opposite sign, as expected from the 
classical Hall effect. MacDonald et al. have studied effects of an injected current on the integer quantum 
Hall state about two decades ago [H(| . They have calculated current and charge distributions and found 
that the charge accumulation occurs in the integer quantum Hall state. The charge accumulation causes 
the energy correction via the Coulomb interaction between the accumulated charges. It is expected that 
the same type of energy corrections exists even in the present highly correlated quantum Hall states. 
However, it has not been studied as far as the present authors know. In what follows, first, we derive 
current and charge distributions in the striped Hall state and the ACDW state using the effective action, 
then, we estimate the current dependence of energy corrections of the two HF states. It is shown that 
the energy of the ACDW state increases faster than that of the striped Hall state as the injected current 
increases. 



4.5.1 Current and charge distributions 

We study current and charge distributions of the striped Hall state and the ACDW state. Let us begin 
with the system with no injected current. We denote the two-point function in the HF theory with no 
injected current as (^(x, t)\l/(x', t))j=o- If a finite electric current is injected, electromagnetic fields and 
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the two-point function deviate from their original values. We define these deviations by 

a(x, t) =A(x, t) - A ex (x), 
Mx,x',t) H^x, - (*t(x,t)*(x',t)) 7=0 , (4.31) 

where a(x, t) and Sp(n, x',t) are unspecified for the moment and will be determined later. The total 
action in the Coulomb gauge V ■ A(x) = is given as 



S to t[A, = J tf^!A 2 (x,t)-l(VxA(x,f)) 2 j 

+ jdtd 3 x^^t) (^- (P + e 2 ^ ,t))2 



*(x,*)<J(z) 



- / did 3 xd 3 a;'* t (x,t)* t (x',t)U(x-x') , I'(x',i)*(x,i)J(z) ( 5(z')- 



(4.32) 



where p, is the magnetic constant, m* is the effective electron mass, and the dot means the time derivative. 
This total action consists of the three-dimensional electromagnetic field term and the two-dimensional 
electron field term. In the Coulomb gauge, the interaction between electric fields is expressed by the 
Coulomb interaction. By applying the HF approximation to the Coulomb interaction part and substi- 
tuting Eq. (|4.3ip into Eq. (|4.32|) . we rewrite the total action as 



(4.33) 



where 



SbmH = J dtd 3 x Qa 2 (x,i) - i-(V x a(x,i)) : 



S*HF[fy,a,1' t ,*] = J dtd 3 x^(x,t) (id t - 



(p + eA cx (x) + ea(x, t)) 2 
2m* 



dtd 3 xd 3 x' 



((¥ (x, t)tf (x, t))i=o + 8p(x, x, t))V(x - x')* V. ^( x '> *) 



- ((* t (x,i)*(x',t)) /=0 +Mx,x',i))U(x-x')* t (x',t)*(x,i) 5(z)(5(z'). (4.34) 

In the expression of Sem [ a ] > the term of the uniform external magnetic field is dropped since it gives only 
the same energy constant to the two HF states. In Eq. (|4.34p . the term including <5p(x,x, t) and the term 
including <5p(x, x', t) are the Hartree term and the Fock term, respectively. As seen in Appendix IB. 1\ the 
Fock term becomes negligible compared to the Hartree term in the long wavelength limit since in the 
momentum space, the Hartree term is proportional to the Coulomb potential V(k), which is 0(l/k), and 
gives a larger contribution than the Fock term for the small momentum k. In the following calculation, 
the deviation of the Fock term is dropped. If we introduce the potential generated by the electron density 
deviation as 



a (x,t) 



,(-e)<Jp(x',x',i) 



47re|x — x' 



S(z'). 



(4.35) 



Shf is rewritten as 



5 H F[o ,a,$U] = J dtd 3 x¥(x,t) (id, 



] t + ea (x,f) 



(p -)- eA cx (x) + ea(x, t))' 
2m* 



dtd 3 xd 3 x' 



(^(x, t)*(x, t)) I=0 V(x - x')* f (x', t)*(x / , t) 



(*t( x , t)*(x', t)) I=0 V(x - i)*(x, t) S(z)S(z'). 



(4.36) 
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This HF action Shf has the same form as the Hamiltonian in the system with the infinitesimal external 
gauge field shown in Eq. (|4.8[) if we apply the HF approximation to the Coulomb interaction part. The 
important difference is that a M (x, i) in the present case represents not the infinitesimal external gauge 
field but the finite gauge field induced by the current flow. Although the meaning of a^(x, i) is different, 
the effective action obtained in the previous section is applicable as long as a M (x, t) is small. 
The partition function is given by 

Z-fttoJrtV,*-******'*. (437) 

Integrating out electron fields and expanding the results up to second order of ao(x, t) and a(x, i), we 
obtain the effective action S c s[aQ,a\ as 

*-jv**«*-»»+«. (4.38) 

The functional derivative of (Sem[o] + So + AS c ff [do, a]) with respect to a(x, t) gives the Maxwell's 
equation for a(x, i), 

( e9 2 _ I V 2 )a(x, t) = (j(x, t)) a 6(z), (4.39) 

where j(x, t) is a current operator and (0(x)) a means an expectation value of an operator 0(x) for 
the system with finite a AI (x, t) (or cquivalcntly finite currents). The solution of this equation gives the 
stationary point of the action with respect to a M (x, t) . We use the action into which the solution of 
Eq. (|4.39p is substituted as the effective action. (j(x, t)) a and Sp(x.,t) = Sp(x, x, t) are calculated from 
the effective action by 

SAS e s[a ,a\ 



<5a(x, t) 
5AS c ff[aQ,a] 
Sa (x.,t) 



--(-e)[p (x)+6p(x,t)]6(z), (4.40) 



where the /Oo(x) is the expectation value of the density operator in the system with no injected current. 
Equations (|4.35|1 . (|4.39|) and (|4.40|) determine ao(x, t) and a(x, t), or equivalently <5p(x, t) and (j(x, t)) a , 
sclf-consistcntly. 

We concentrate on the finite system with the static injected current flowing in the y-direction and 
depending only on x. The lengths of the 2D electron system in the x-direction and the y-direction are L x 
and L y , respectively (Fig. 14. 6p . In this case, the electron density also depends only on x and Eqs. (14.351) 
and (|4.39|) give the following solutions at z = 0: 



1 



2ire 



L x /2 



ao(x) = — - — / dx' In \x — x'\(— e)Sp(x'), 



-W2 



a y (x)=^- ( ^ dx'ln\x-x'\(j y (x')) a . (4.41) 



2tt 



-Lx/2 



As shown in the previous section, the effective action can be divided into the non-perturbed part and the 
correction part due to currents, and in the long wavelength limit, the correction part AS'off is given by 

dx(-e)p a a (x) dx a^(x)K£ v (d x )a v (x), (4.42) 

where po is a uniform part of the density, and T is the total time. Kq V (8 x ) is the Fourier transformed 
form of the response function obtained in the previous section. Substitution of Eq. (|4.42[) into Eq. (14.401) 
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Figure 4.6: Schematic view of the 2D electron system in a magnetic field with the injected current. The 
current flows in the y-direction and has only the cc-coordinate dependence. 



gives 8p{x) and (j y (x)) a as 



(-e)Sp(x) =K°°(d x )ao(x) - K° y (d x )a y (x), 
( Jy (x)) a =K y0 (d x )a (x) - K yy (d x )a y (x). 



(4.43) 



Equations (14.41j) and (|4.43p determine the current and charge density distributions up to an overall 
constant. The overall constant is determined by requirement of the following constraints: 



W 2 



dxjy (x) 



-L x /2 



ix/2 



L x /1 



dxSp(x) = 0, 



(4.44) 
we obtain 



where J is a total current. Using the explicit form of the response functions derived in Sec 
the integral equations to determine the current and charge distribution. 

The same type of the integral equations has already been solved for the integer quantum Hall state 
50l . 51, 52, 53| . Their results are summarized as follows: 



1. In Eq. (|4.43p . the terms including the vector potential a y {x) give a very small effect in the integral 
equations compared to the terms including the scalar potential ao(x) and the vector potential terms 
are negligible in a good approximation. 

2. The analytical solution of the integral equation without the vector potential term is obtained by 
means of the Wiener-Hopf technique. 

3. ao(x) = const x In \(x — L x /2)/(x + L x /2)\ is the good approximate form of the analytical solution 
except near the edge and the constant coefficient is determined from the constraint for the total 
current. 

The same results hold in the present case. The integral equation for the potential is given by 

,L x /2 



ao(x) = — 7 / dx' In \x — x / \d x ,ao(x'), 

J-L x /2 



(4.45) 



where 7 = (1 + ft\pB j 'v)a x V y /2-k€lo c {ft = for the striped Hall state). 7 has the dimension of length 
and is very small for the magnetic fields of the order of several tesla in the quantum Hall regime. For 
example, if e = 13eo, m* = 0.067m e (these are parameters in GaAs), and ft = 0, then 7 is of the order of 
10 -8 m. The current and charge distributions are obtained from ao(x) by 



(-e)Sp(x) = 2neyd^.a (x), (j y (x)) a = -a^d x ao(x). 



(4.46) 
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Figure 4.7: Potential distribution ao(x). The first derivative of ao(x) gives the current distribution and 
the second derivative gives the charge distribution. 



The approximate solution of Eq. (|4.45|) is given by (Fig. I4.T[) 

for |x|<^-7, (4.47) 

with a linear extrapolation of ciq to ±IRh/2 in the interval within 7 from the edge, where a = 
IRh/2(1 + ln(i x /7)) and Rh = l/(Jxy is the Hall resistivity. One may verify that Eq. (|4.47l) is in- 
deed the approximate solution of the integral Eq. ()4.45j) . by substituting Eq. (|4.47j) into Eq. (|4.45|) and 
performing one partial integration. 



ao(x) = a In 



x - L x /2 



L x /2 



4.5.2 Energy corrections 

The energy correction due to the injected current per unit space-time volume is calculated from the 
effective action by (5em[i] + AS c ft [a>o, a\)/TL x L y . Since in the present case of v* = 1/2, the area 
occupied by one particle at the uppermost partially-filled LL is 2a 2 (here, the vNL constant a is written 
explicitly), the energy correction per particle 5E is given by [(SemM + AS c ff[ao, a])/TL x L y ] x 2a 2 . 
Substituting Eqs. (|4.41[) and (|4.43[) into this expression, we obtain the energy correction per particle 
given by 

P 2 pL x /2 

5E[I] = -- — — / dxdx'5p(x) In \x - x'\5p{x'). (4.48) 
2ireL x J~l t /2 

Substituting Eq. (|4.46|) into Eq. (|4.48j) and using Eq. (|4.45() . we rewrite the energy correction as 

5E[I] = 1 xy) dxao(x)d 2 x a (x). (4.49) 

Substituting Eq. (14.4711 into Eq. ()4.49p and performing the x integral, we obtain the final result; 

SE [I] = x ^§ L x l\ (4.50) 

(ln(2/&) + l) 2 V j 

where b is a dimensionlcss constant given by 6 = 7/(£ x /2)(<C 1). This expression depends on the filling 
factor, the magnetic field strength and experimental parameters. Since the actual filling factor includes 
the spin degree of freedom, we use i/ ex = 2lo + v* for lower spin bands and v ox = (21q + 1) + v* for upper 
spin bands instead of v. The magnetic field strength is related to the filling factor by B = hn e /ev cx (n e 
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v cx stripe parallel perpendicular I c 
5/2 325.0 330.6 335.6 0.041 
7/2 204.4 206.6 209.0 0.065 
9/2 144.7 146.1 147.6 0.040 
11/2 109.9 110.7 111.6 0.053 
13/2 87.44 88.08 88.68 0.047 



Table 4.3: Values of the coefficient C in units of nA -2 and the critical current Ic in units of nA. parallel 
is the case in which the stripe direction is parallel to the current, perpendicular is the case in which the 
stripe direction is perpendicular to the current. 




Figure 4.8: Energy differences AE[I] between the striped Hall state and the ACDW state. The results 
at v = 5/2, 9/2, 13/2 are shown. When Ae[I] is positive, the striped Hall state has lower energy. 



is an electron density). For example, if n e = 2.67 x 10 15 m 2 , then the magnetic field strengths are 4.42 
T (i/ ex = 5/2), 3.15 T (zy cx = 7/2), 2.45 T (zy cx = 9/2), 2.01 T (v cyi = 11/2), 1.70 T (v cx = 13/2) and so 
on. We use e = 13eo, m* = 0.067m e , n e = 2.67 x 10 15 m, and i x = 5x 10~ 3 m in order to estimate the 
values of energy corrections, which are the parameters used in the experiment by Lilly et al. 0. Then 
the energy correction is given by SE[I] = C x I 2 (q 2 /Ib) with the coefficient C shown in Table l4~3l 

As shown in Sec. 14.31 in the system with no injected current the energy of the ACDW state is slightly 
lower than that of the striped Hall state. The differences of energy per particle AEq are 9.3 x 10~ 3 
(la = 1), 2.3 x 10~ 3 (lo = 2), 1.4 x 10~ 3 (Iq = 3) and so on in units of o 2 /Ib- When the finite current is 
injected, charges are accumulated in both edges with the opposite sign. The accumulated charges give the 
energy corrections SE[I] which depend on the value of current I. Including these corrections, the energy 
difference between the striped Hall state and the ACDW state AE[I] = — AE + (5Eagdw[I]~ ^stripe [-f]) 
varies depending on I. The current dependence of AE\I\ is shown in Fig. 14.81 In Fig. 14.81 only the parallel 
case is plotted for the ACDW states since it has a weaker current dependence than the perpendicular 
case does. The signs of the energy differences change at the critical values of current I c . The critical 
values are shown in Table EOfl The critical values are about 0.04 - 0.05 nA for I — 2 and I = 3 LLs. The 
current used in the experiments 0, [2] is above 1 nA and is much larger than the critical value. At 1 
nA, the energy differences Ae[I] become 1.4 (y = 9/2) and 0.63 (v = 13/2) in units of q 2 /Ib which are 
much larger than the original energy differences at a zero injected current. Hence, the striped Hall state 
becomes the lower energy state and should be realized in the experiments. We note that for / = 0.1 
nA, the density deviation around edges is of the order of 0.001% of the mean value of the density, which 
ensures that the electromagnetic fields induced by currents are kept small enough for the current of the 
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order of 0.1 nA. 

4.6 Summary 

In this chapter, we have investigated the effect of the finite electric current on the striped Hall state and 
the ACDW state in the system with no impurities and no metallic contacts within the self-consistent HF 
approximation. First, at zero injected current, the energies of the two HF states have been calculated 
and we have confirmed that the ACDW state has slightly lower energy than the striped Hall state 
as pointed out in Ref. (28l |. Next, within the same HF approximation, energy corrections due to finite 
injected currents flowing in the stripe direction have been calculated and the energies of the two HF states 
including the energy corrections have been compared. For this purpose, we calculated electromagnetic 
response functions, and using them, current and charge distributions were determined for both states in 
the system with the injected current. It has been found that the charge accumulation occurs around both 
edges with the opposite sign, just as in the case of the integer quantum Hall state studied by MacDonald 
et.al and other authors [H(J, El, HH . We hope that current and charge distributions will be observed 
in experiments for anisotropic states. The charge accumulation results in the energy enhancement via 
the Coulomb interaction between the accumulated charges. The energy enhancement was estimated from 
the current and charge distributions. It is found that the energy of the ACDW state increases faster 
than that of the striped Hall state does as the injected current increases. Hence, the striped Hall state 
becomes the lower energy state when the current exceeds the critical value. The critical value is estimated 
at about 0.04 - 0.05 nA. The current used in the experiments for the anisotropic states @, H3] is above 
1 nA. This result suggests that the striped Hall state is realized in the experiments. In addition, the 
striped Hall state has the anisotropic Fermi surface, which naively explains the experimental features of 
the anisotropic states, i.e., the anisotropic longitudinal resistivities and the un-quantized Hall resistivities. 
Hence, we conclude that the striped Hall state is realized in the experiment rather than the ACDW state 
and predict that the ACDW state is realized if an experiment is performed with a current smaller than 
the critical value. 



52 CHAPTER 4. EFFECTS OF AN INJECTED CURRENT ON HIGHLY ANISOTROPIC STATES 



Appendix A 

von Neumann lattice formalism 



A.l Wave function of the von Neumann lattice basis 



By using the von Neumann lattice basis introduced in Chapter [51 the Hilbert space of a one-particle state 
of a free electron in a magnetic field B = (0, 0, B) is spanned by \l, N) = |/;) <E> \a m . n ). The wave function 
in the coordinate space i/';,n(x) = (x|Z,N) is given by 



^,n(x) = e 
where </>z,m(x) = (x|i, M) with 

\l,M) 



Z7r(m+n+mn) g — 7r| 2 m T1 1 /2a 



2 /9n 2 \— * 



1 / \ZnZ m n 



M 



^,m(x), 



(t \ \l crt \M 

1 A) [ B) |0,0), L A |0,0) = L B |0,0) =0. 



The explicit form of (^a/ (x) is obtained by using the coordinate expression of La and Lb-, 



La 



V2h 



-(x — iy) — il 



(- 


-«— Y 


\dx 


dy] 



y/2h 



(-- 


■ 9 M 


\dx 


^1 



(A.l) 



(A.2) 



(A.3) 



These expressions are rewritten by using the complex coordinates z = (x — iy)/ls and z* = (x + iy)/ls 



as 



9 



dz 
c 

dz 



/4 (z * _ 2 |. )e l^ 

9 ui 2 //i 1 _ui 2 //i, „ d 



Lb = V2e-^/^e^/\ L B = ^^(z - 2^)e^l 

V2 dz* 



(A.4) 



where the derivations act on not only the exponential function but also the wave function on which these 
operators act. Using these expressions and performing some calculations, we obtain the explicit form of 
^.a/(x); 

where C;.m is a normalization factor, and L^(x) = (e x x~ n /m\)(d/dx) m (e~ x x m+n ) is a Laguerre poly- 
nomial. In Eq. (|A.5|) . the polar coordinate representation of z, z = r(cos9 — ismO)/lB, is used. If I = 0, 
0o. m (x) is given by 

0o,m(x) = -==!=- 



V2w2 M M\ I b V ^ 
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Equation (|A.6p gives 

0ON(x) = — e in ( m + n + mn ) e -i(ir/a)(mr 3 y-nx/r 3 ) e -(Tr/2a 2 ){(x-amr 3 ) 2 + (y-an/r a ) 2 } (A 7) 



|<fo n(x)| 2 = -L e -^/^H(—mr s y + (y-an/r a y}^ ^ 

1 * 



1 

Hence, for Z = 0, the probability density is localized at the lattice points, a(mr s , n/r s ). A similar argument 
can be applied for the case of I ^ 0. 



A. 2 Explicit form of the inner product (a p \a p >) 

By using the momentum representation of |a m , n ), \ce p ) = J2 m n e mP:L+mPy \a m ,n), the inner product is 
given by 

(Op|<V)= Yl e " JP ' N+iP '' N '(«m,nl«™',«')- (A.9) 



m.n.rn' .n' 



Since the {(X m n \a m t n *) is a function of (m — m!) and (n — n'), Eq. (|A.9[) can be rewritten as 

(a p |cv) = a(p)(27r) ^ S 2 (p p' 2ttN), 

N 

a(p) = J2 e~ v * M - tPyN (uM,NM, (A.10) 

where the relation, Y^=-oo eWn = (2 71 ") X)^-oo ^(p — 27m), is used. In what follows, we derive the 
explicit form of a(p). 

From Eq. (|2.15|) . (a m ,„|ao,o) is given by 

TT 

(am,n|ao.o) = exp(-—x Ax + «7r(m + n)), (A. 11) 

where 

■-(:)■ (A - l2) 

By using the identity, 

1 = -(detA)" 1 / 2 / d 2 s exp(--(s-aa;) t A- 1 (s-aa;)), (A.13) 

with 

.-(*), detA = 2, A- = i(^ '), (A.14) 

the index of the exponential function in Eq. (|A.13[) can be transformed from the quadratic function to 
the linear function in m and n as follows, 

1 f°° TT TT 

(om,n|ouo) = — 7= / d 2 s exp(— — (s — iAx) 1 A^ 1 (s — iAx) ) exp(— — x* Ax + Z7r(m + nj) 
2\/2 J-oo 2 2 

d 2 s exp(— — ■s t A~ 1 s + iTT{ms x + ns y ) + i7r(m + ri)). (A. 15) 



2^2 7-00 2 
Substitution of Eq. ([A~T5)) into Eq. ([ATO]) yields 



Q:(p) = / e -(ir/2)s t A- 1 sT^ e -ip :c m-ip y n+iir(sx+l)m+iTr(s y + l)n^ (A. 16) 

2^2 J oo m ^ n 
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The summation of m and n gives the delta function. The integration of s x and s y rewrites a(p) as 

P x /tt + 2m — 1 



(p) =\/2 J^exp 

m , n 

=V2 exp 



_J(^ + 2m-l,^+2n-l)A- 1 . 

2 TT 7T V Py/ 71 + 2n — 1 



7T f 1 /J3; 



4 I 7~ V 7T 



2 (Py 



T ^ --1 ^-1 -^ --1^-1 



2 V TT 



Py 



x ^ exp 

n I 

x ^exp 



7T 


M- 


i 

^ s 


rpx_ 

V TT 


- 1 


7T 


-H 


-ir 2 s 


<Py_ 

V 7T 





(A.17) 



By using the relation between the summation of an integer n and the definition of a Jacobi's theta 
function of the first kind 6\{v\r), 



e iir(l-T+2v)n+iTTTn 2 _ _ ip"™^ /4-v) ^ ^ ^ 

n— — oo 

the summation of m in Eq. (|A.17[) is rewritten as 

-ie-^'^OMn), with Wl = * (p^+fr^), ri 
and the summation of n in Eq. (|A.17|l is rewritten as 



ie-^/^eMn), with «a=-^-(R B -irJp w ), r 2 



,•2 ■ 



In addition, by using the relations between theta functions, 

0i(v|t)=-0I(-v*|t), 
^( w |r) = e^V^/V^^^/Tl - 1/r), 

here we provide that r is a pure imaginary, we obtain the explicit form of a(p)) given by 

a(p)=/3(p)/3*(p), 

where 

/3(p) = (V^r,) 1 ^-^) 2 /^ I P* + ir2 sPv 

\ 2tt 

A. 3 Density operator on the von Neumann lattice basis 

On the von Neumann lattice basis, an electron density operator p(k) = J d 2 x W (x)^f(x)e x is expressed 
as 

d 2 p d 2 p' 



(A.18) 

(A.19) 

(A.20) 

(A.21) 
(A.22) 

(A.23) 
(A.24) 



'BZ 



(2ttY 



p<m = E / 

The integral over x is performed as 

J d 2 x(/,p|x)(x|/',p)e lk - x = (Z,p|e tk '( x+ «) / d 2 :. 



tf(p)Mp') J dM*, P |x)(x|Z',pV k - 



(A.25) 



\l',P , ) = (.fi\e^\fi')((3 P \e^\(3 pl ), (A.26) 



• r j(i 
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where the operation form of x, x = X + £, and the identity 1 = Jc? 2 a;|x)(x| are used. In Eq. (|A.26j) . 
(/j|e k '£|/;<) = /°;/(k) is given in Appendix IA.41 Whereas, (/3 p |e 4kx |/3 p ') = F(p,p',k) is calculated in 
what follows. 

By using the relation |/3 p ) = |a p )//3(p), F(p, p',k) is expressed as 



F(p,P , ,k) = 



1 



(a p \e^ x \a p ,) 



(A.27) 



Using the relations X = (a/2y / 7r)(is + L B ) and Y = — i(a/2y / 7r)(Ls — L B ), we rewrite e 4k x as 

e ik-X _ e (a/2y^)(ik T +k v )L B + (a/2^)(ik :c -k y )L i B _ e ~{a 2 /?,Tr)k 2 e (a/2 x /¥)(ifc x -fc !/ )L B e (a/2 v /ff)(ifc a .+fc„)L B rp^ 28) 

here, the Cambell-Hausdorff formula Eq. (|2.21|) is used. By using Eq. (|A.28|) and the relations, 

^s|^m,n) — ^m,7i|^m,7i)' ( ^771,71!-^ (^m,n | ~ ^m,n> (A. 29) 

the matrix element (a p |e lk ' x |o! p <) is calculated as 



/ 1 ik-Xi 

(a-a e a r 



E 



-(l/87r)/c 2 £ ,(l/2)(ife x -fc H )z^ „ + (l/2)(ifc x +fc B )z m , „, -ip-N+ip'N' 



N,N' 



p -(l/87v)k 2 ^ e -i{(l/2)(px-(rj2)(kx+ik y ))} m -i{p y + (l/2r s )(ik x -k y )}7i( 
N 

x y^ e HP»+(^/2)(fc»-ifc a )}m , +i{ P ; + (l/2r s )(ifc a ,+fc y )}n'| Q ,^ ) ^ 



N' 

- e (l/8ir)* (a px _( r3 /2)(fe x +ifc H ),p M +(l/2r s )(ifc x -fe H )|"p^ + (rV2)(fex-ifcy),P^ + (l/2r a )(ifc x +fe s )) 



3 "/bx - y (** + iky), Py + ^-(ik x - k y ))(2n) 2 S 2 ( P - p' - k - 2ttN), 

rs N 

(A.30) 

where we set a = 1 and define /(p) = |/3(p)| 2 . Since the explicit form of /(p) is given by 



/(p) = (V2r s )e-^ 2 ^ 



Px - »r 2 Pj; 
2tt 



Px + 

2tt 



(A.31) 



the function f(p x — {f s /2){k x + ik y ), p y + (l/2r s )(ik x — k y )) in Eq. (j A.30|) is given by 
f(Px ~ -f(k x + iky), Py + ^{ikx - k y )) 

= e fc2 /87r e -(r 3 p B ) 2 /47r e -^(p ! ,-fc B /r 3 ) 2 /47r e -( 4 r s /47r)fc x (2p ! ,-fc ! ,/r ii ) 



Px ~ iriPy 
2ir 



(p x -r s k x ) + ir 2 s {p y - k y /r s ) 
2ir 



(A.32) 



Hence, 



(a p |e ik - x |ov) = (3*{p)(3(p - k) e -(^/ 4 -)M2p H -fc y /r 3 ) (27r)2 S 2 (p -p'-k- 2ttN), (A.33) 



N 



and 



(/?p|e lk ' X |/V) = e -(^s/irr)kA2p y -k y /r a )^2 § 2 (p p' 2ttN> 



5 #(P',N) 



(A.34) 



N 



where the boundary condition for /3(p) is used. As a result, the density operator p(k) is expressed on the 
vNL basis as 



^)=E/ BZ ^ t cp)^(p , -^« 



(k) e~ ' irs / 47r ) fc * ( 2 Pb ~ k v / r * ) _ 



(A.35) 
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A. 4 Landau level matrix elements 



The matrix elements (Zi | e zq "^ | Z2) are given as follows 

h-i 



{h\e^\h) = 



ft e-^;^ (£) for h > h 

for h = k, 



(A.36) 



where L\ (x) is a Laguerre polynomial. For I 1 = 0, the Laguerre polynomial is given by 





n=0 


J! 

(1 — n)\n\n 


L° (x) 


= 1, 




L\{x) 


= 1 — X, 




L° 2 (x) 


= 1 - 2x + 


1 2 

2*' 


L° 3 (x) 


= 1 - 3a; + 


3 9 1 
-x^ - -x 6 . 
2 6 



From Eq. (|A.36|) . the matrix elements ; (q) defined by Eq. (|4.10|) are given by f? l (q) = (l\ 
f* (Ja (q) = iu;<A v <Zi|e^%) and /f^q)* - -^c^ile^). Note that {/£^(-k)}* = 
holds following from its definition. 

The values of j 2 (0) and its derivatives are given by 



(A.37) 
/f 2>il (k) 



<,z a (?) 



,Jh±±8 h<h+1 fovl 2 >h 



dq y 



9 =0 



for li > I2 
for /1 = 1-2 



A/ l -rr s h,h+i br l 2 >h 



q=0 



for li > h 
for /1 = Z2 



(A.38) 



(A.39) 



dq v 



dq x 



-^ ^1 + 1)0 1+2)^^+2 forfc>ii 

-i^^h{h - l)8 hM+2 for Zj > Z 2 

-*&Gl + 5) for/ 1 = / 2 

^y /(li + l)(I iT2)^ 1+2 for/ 2 >^ 1 

-i^^h{h - l)8 hM+2 for Ji > h 

I + |) for h = h 



A. 5 Hartree-Fock Hamiltonian 



(A.40) 
(A.41) 



In the HF theory, the Coulomb potential term given by Eq. (|2.35p is approximated by Eq. (|2.36j) . In 
the summation in Eq. (|2.36| . both the Hartree term and the Fock term are proportional to the density 
operator p; 3 .i 4 (k). This is a characteristic of the quantum Hall system 54], 55|. In what follows, we derive 
Eq. (|2"3o) from Eq. (2J55) . 
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A. 5.1 Hartree term 

The Hartree term is given by 



V Hartr e e = £ I 4% I ^( k )/l 1 ,Ia( k )/«°„I 4 (- k )<rt I .»»( k )>W3,u(-k)- (A.42) 
Transforming the variable k to — k and using the relation V{— k) — U(k), we obtain 

V Har t r cc= £ / TjSa / ^( k )/h,i a (- k )i?a,u( k )<P'i,fa(- k ))PJ3,u( k )- (A.43) 



l\ ,^2 



A. 5. 2 Fock term 

Fock term is given by 



VFock = - £ SwfL ^y^<^(p)Mp'H(p' - k)6 i4 ( P - k)i££ i{s>u (k), (A.44) 



with 

<t W J k ) = ^(^/^^(-^^^(^e^^)^^^). (A.45) 

For simplicity, we define V/ by VFock = — z 2 i 3 z 4 Vl an d calculate V/. 
If we insert the identity 

1 = y ^^(27r) 2 <S(q - ^e^ 2 ^^^"^'-^ 2 ^^^-^) (A.46) 

into V/, we obtain 

v " 1 1? / (0 1 W <4 <* < p > t "< p '» !>, < p ' - k > 4 < p - k «,w.M 

x e -(V 27r )<?^(PB-pS,)~( i / 2,r )'J!'(Px-Px)^27 r ) 2 ^ 2 (q _ ] t ) e ( i / 2T ) fc y(pL-p^) > (A. 47) 



Dividing the integral range of the k integral by 

72 



f d z k , v-^ f d k 

iM flk, =?I,M f(k ^' Nl 



(A.48) 



and using the relation derived from the boundary condition Eq. (|2.29| 

& f (p' k + 2ttN)6( P - k + 2ttN) = 6 f (p' - k)6(p - k)e i(p --^ )JV « , (A.49) 

we rewrite V; as 

Vi = L d2 X^ V (& ' i(p)6 ' 2(p,))& ' 3(p ' ■ k)Mp ■ k x^ i3) u(q) 



x e -(»/27r) 3<e (p v -p;)-(i/27r) 9v (p;-p ai ) ^(27r) 2 c5 2 (q k + 27rN)e (4/27r)fc! ' (p -~ Pa:) . (A.50) 

N 

Using the relation 

Y J {2n) 2 5 2 (q-k + 2nN) = ^2e-'^- k ^ N y+^ q y~ k ^ N ^ (A.51) 



N N 
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and rewriting the variable k as p' and p' as k, we obtain 



Vl ^lwr /B Z ^ (6|i(p)i,2(k)> ^ (k " p/)b(p " p>& ^® 

x e -(i/^)Qa,{p v -ky)-(i/2^)q v (k x -pa S )+(i/2^)p' y (k x -p x ) e -i(q x -p' x )N v -H(q y -p' y )N x ^ (A. 52) 

Transforming k to k + p and p' to — p', we obtain 

v '=E/(0 J bz rf2fc ff V ( b t Crib, (p - k M 3 (-k + p + P')bu (p + p'K^.u (q) 

x e -(«/ 27r )9xfcB + ( i / 27r )9B' i: x e (i/27r)p^fcx e -i(gx+Px) Ar B +P„)^x (A. 53) 

Using the relation 

6, 2 (p - k)fej 3 (-k + p + p') = b h (p - k - 2^N)6[ 3 (-k - 2ttN + p + p')e ip « , (A.54) 
we rewrite Eq. (j A.53|) as 

V <=E / BZ (02 / (0^^ 1 (p)^(p-k-2 7 rN ) ) & t 3 (- k -2 7 rN + P + p')^(p + p') 

x yfock^ ^ ^ e -{i/2K)q x (k y +2TTN y ) + {i/2K)q y (kx+2KN x ) e {i/2T!)p' y {k :c +2T!N :c ) (A. 55) 

Using the relation 

?I(0 F(k+2 - N) = /(§^ (k) ' (A - 56) 

and transforming p' to p' p, we finally obtain 

^^^(q)^™^^- 9 ^^- (A.57) 
Thus, the combination of the Hartee term and the Fock term gives 

V hf= E /^<l, i3 , i4 ( k )(^^(- k )>^,u(k), (A.58) 

with the HF potential 

< F l2 , 3 , h (k) = ^(k)/0 )i2 (-k)/f 3iU (k) - | ^^(kO/^^-kO/^^kOe-^/ 2 -^^^-^^). (A.59) 
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Appendix B 

Highly anisotropic states 



B.l Deviation of the Hartree-Fock Hamiltonian due to currents 

In Sec. 14.51 the deviations of the magnetic field and the two-point function arc taken into account in the 
long wavelength limit in order to consider the current effect on the HF states. The deviation of the two- 
point function is caused by the deviation of (p; 1) ; 2 (— k)). We only consider the deviation at the partially 
filled LL Iq since it would give the largest contribution in our calculation and denote it as Spi (— k). Then, 
the deviation of 7Yhf is given by 

Whf = £ J pCwA^A,^). (B.l) 

In the long wavelength limit, 5pi {— k) is relevant only for the small momentum. When we expand 
v h F i 2 1 U ( ) w ith respect to k, the largest contribution comes from the lowest order term with respect to 
k. For each set of the LLs (I3, 14), the Hartree term of the HF potential gives the lower order term with 
respect to k since the Hartree term has V(k) which is 0(1/ k). Hence, the Hartree term gives the main 
contribution and the Fock term is negligible in the long wavelength limit. 



B.2 Calculation of K 00 for the ACDW state 

When po = Py = an d p x — » 0, the response function K^ ) (p x ) is Taylor expanded with respect to p x as 



K° °( Px ) = K o (0)+p x d p „K° o (0) + P -^dlK™{0) + 



(B.2) 



The first and second terms become zero. The third term includes the corrections from the inter-LL term 
and the intra-LL term. The inter-LL term gives the same expression for Kq° as that in the striped Hall 
state. The intra-LL term gives the extra correction AKq°(p x ) given by 

AK°°( Px ) =KXp% (B.3) 

where k is given by 



d 2 p' 



1 



2 J RBZ (2^2e(p') 9p 



A{p + p')A(p') + Rc(B(p + p')B*{p')e-^ 2 ^ 
e(p + p')e(p') 



(B.4) 



Px=0 



/3 in Eq. (|4.30|) is defined by (3 = — k x vlo c /(JxyVB. The finite k is the result of the band formation 
at the partial filled LL, while the band structure is generated by the density modulation of the ACDW 
state in both directions. Hence, it may be considered that the finite k reflects the remaining density 
modulation effect of the ACDW state in the long wavelength limit. 



Gl 
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